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7.1  Introduction 

Using  altimeter  data  for  quantitative  information  about 
sea  ice  is  an  idea  with  relatively  few  proponents.  Often  the 
poor  spatial  sampling  given  by  the  altimeter’s  single-point 
measurements  along  widely  spaced  ground  tracks  and  the 
difficulty  of  interpreting  altimeter  pulse  echoes  over  ice  are 
cited  as  reasons  to  avoid  altimetry.  It  is  becoming  evident, 
however,  that  altimetry  may  be  able  to  make  unique  mea¬ 
surements.  For  instance,  the  altimeter  data  record  reveals 
the  presence  of  small  areas  of  open  water  within  the  pack  at 
concentrations  too  low  to  be  detected  by  a  passive  micro- 
wave  sensor  such  as  a  Special  Sensor  Microwave/Imager 
(SSM/I)  and  too  small  to  be  resolved  by  satellite  synthetic 
aperture  rada:  (SAR). 

At  high  latitudes  the  net  laid  down  by  a  satellite 
altimeter’s  ground  track  becomes  dense,  and  spatial  sam¬ 
pling  (over  a  period  of  days)  can  approach  the  resolution  of 
SSM/I.  A  satellite  altimeter’s  along-track  resolution  of  less 
than  a  kilometer  is  better  than  that  of  SSM/I,  and  the 
altimeter’s  sensitivity  in  detecting  the  ice  edge  is  greater. 
The  low  data  rate  makes  it  possible  to  process  altimetry  in 
near-real  time.  For  these  reasons,  an  ice  edge  product 
derived  from  Geosat  altimeter  data  has  been  used  opera¬ 
tionally  at  the  Navy/National  Oceanic  and  Atmospheric 
Administration  (NOAA)  Joint  Ice  Center  to  supplement 
infrared  imagery  for  estimating  ice  edge  position  [Hawkins 
and  Lybanon,  1989J. 
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Objections  concerning  the  difficulty  of  interpreting  al¬ 
timetry  over  ice  are  more  difficult  to  ireet.  in  part  becau.se 
of  the  paucity  of  ground  truth  or  imagery  coincident  with 
altimetry.  In  this  chapter,  we  cover  some  of  the  rc.search 
into  understanding  altimeter  backscatter  over  sea  ice  and 
interpreting  altimeter  data  in  terms  of  parameters  such  a.s 
ice  type,  concentration,  edge  location,  and  wave  penetra¬ 
tion  in  the  marginal  ice  zone. 

7.2  Instrument  De.sckii>tio,n 
7.2.1  Instrument  Design 

Satellite  altimeters  have  evolved  steadily  from  Skylab, 
which  in  1973  demonstrated  the  concept  of  altimetry  from 
space.  Development  continued  through  GEOS-3  in  197.5 
until  the  1978  launch  of  Seasat,  which  carried  the  first 
altimeter  capable  of  measuring  sea  surface  topography 
with  an  error  of  only  a  few  centimeters.  Precision  i  n  range 
or  height  measurement  is  obtained  through  a  design  that 
takes  advantage  of  the  ocean’s  near-Gaussian  w'ave  height 
distribution  and  uniform  dielectric  characteristics.  Over 
ice,  where  surface  height  distributions  and  dielectric  char¬ 
acteristics  vary  sharply  from  place  to  place,  measurements 
of  height  and  backscatter  are  often  misleading.  Therefore, 
it  is  necessary  to  understand  how  altimeters  arc  designed 
to  operate  before  one  interprets  altimeter  backscatter  from 
sea  ice.  Altimeters  planned  for  the  1990’s  will  refine  rather 
than  change  the  Seasat  design,  upon  which  the  following 
description  is  based.  Hereafter,  similar  altimeters  such  as 
Geosat  and  ERS-1  will  be  called  ElRS-1  class  altimeters. 
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7.2.2  Fuls,'  l.t/irJfil  .Altifut-t!  ': 

Satellite  aUinu'ters  u.-,e  ;>ul,Tf  s',  u.llh-lit!uliAlwe«Miselr> .  A 
wide-boam  (with  a  :i-dB  halt-asidth  o!  Ot!'  s'lir  St‘a:.ai' 
relatively  loiiit  uin  the  urder  ot  inicm.-iecoiRl.s^  pulse  is 
enutled,  and  the  pulse  echo  i.>  [irocessed  m  a  svay  etpiivak-nt 
to  rneasurine  the  travel  time  to  tlie  suit'ace  and  hack  of  a 
much  shorter  pulse  only  a  tew  nanoc^ccursde  li-my  'I'his 
technique  ot’pulse  coni  pro, -vS  ion  stives  hiqh  range  resolution. 
I  For  a  full  espLin.ition  of  pulsi-  cornpresaio!)  on  -.atcJhle 
altimeters,  see  Chelton  et  al.  il!),k;)!.K  The  area  on  the 
surface  illuminated  by  the  pulse  grows  with  time  until  the 
trailing  edge  of  the  pulse  leaves  the  lowest  renecting  (Xiinl.s 
at  nadir  I  Figure  7-!  K  The  tnaximuin  area  simultaneously 
illuminated,  then,  depends  on  the  surface  height  distribu¬ 
tion  and  etl’ective  icumpressedi  pulse  length.  'I'his  area  is 
usually  referred  to  us  the  pulse-iimitedfootpnnt<PLFi.  Tor 
Seasaf .  thediameter  jfthe  PLT’",-i- t  dlo'T'oi  a.  •'■oothsea 
and  grew  to  7.7  km  for  a  significant  waive  height  *S\VH )  of 
10  rn.  As  the  pulse  continues  to  move  out  over  the  surface, 
the  footprint  becomes  an  annulus  of  constant  area  and 
increa.sing  circumference.  It  i,s  important  to  note  that  due 
to  irregularities  in  surface  roughness  and  slope,  the  true 
illuminated  area  at  any  time  may  be  inx'gular  in  shape  or 
discontinuous.  Over  ridged  sea  ice,  for  instance,  an  otf- 
nadir  ridge  sail  may  be  high  enough  to  be  illuminated  at  the 
same  time  as  the  ice  near  nadir. 

Altimeters  that  use  beam  width-limited  geometry  emit 
a  narrow  beam  and  do  not  rely  on  pulse-compres.sion  tcch- 
njques  for  accurate  range  measurements  Footprint  size 
depends  only  on  beam  width  and  satellite  altitude.  Thi.s 
geometry  is  not  used  for  satellite  altimeters  because  precise 
antenna  pointing  capability  is  a  necessity  and  the  antenna 
size  and  gain  required  to  generate  a  narrow  beam  are  large. 
However,  use  of  a  beam  width-limited  altimeter  would 
simplify  interpretation  ofsignals  over  ice,  because  footprint 
area  remains  constant  rcgardle.ss  of  the  surface  height 
distribution  (Rapley  ct  al.,  19831. 

(aj 


Pulse 


Time 

Fig.  7-1.  (a)  Area  illuminated  by  a  pulse  bitting  the  smooth  ocean’.s 
surface  and  (bi  the  growth  of  the  footprint  .area  with  time  for 
smooth  and  rough  seas.  The  PLF  is  largerfor  meter-scale  roughness, 
since  the  time  required  for  the  nulse  to  travel  from  crest  to  'rough 
oxcoods  the  off.  cliv.  puT.,  'ongth.  Graph  is  not  to  scale. 


(  to  tiu.liii  the  --.Il.'llite  t!u  I  I  !u.>  et  lli'  .  iiiitt.  d  jHlise  i^ 
ditfvrenced  with  a  ili-i  .iinpuic  -'.KKutl  i  lie  !i*  qui’ro.  \  el  the 
emitted  pulM,-  chalices  i.Ae!'  ttiv  pul  e  ieisyth  iii  a  hiiear 
chirp  III  I  amp  I'he  tol.il  cb.m; a  so  i  n  (jueiu  >,i\  t-r  !  in- 1  arn  ji 
I.S  the  baiulw  idl h uf  ih.e  i;s.--tniiisi-ii!  I'lie  del  .u'.ipiiir  >>;.-.nal 
'su-Ciilled  hiA.tUM'  till*  elleci  el  thl-  priHi  -Mii;’  ^tep  r-  to 
remove  the  ramp  ui  liequeiicy  inun  thi.-  piil.-.e  i-iliu  i- 
gellei  aled  at  .i  delas  tUoe  /..  deti-rm ilh,  d  iw  the  onijiKird 
tracking  uiul  Tla.'  time viirre- jiund.-  In  the expi-i  ted  .irnvai 
time  at  the  recei\er  ef  the  reflected  pu!m'  froni  Hican  .-ea 
level.  The  re'-iiltmg  .'ignal  i.-  a  mixture  id  return.-  from 
individual  ponil  .--catti.  rer.-.  Kchi»e-  iiotn  w  .necre.-t-  arni  e 
fir.-t  and  are  therelore  .shifted  earlier  in  time  IVum  tracking 
time/,,.  For  a  pul.-e-hmited  altimeter,  returns  -liil'ted  later 
than  /,,conu‘  not  only  from  jioint.-  lu  low  mean  sea  level,  liut 
also,  as  the  fiKitprint  expand.-,  from  .-catlerer.- distant  from 
nadir  ifthe  height  di.strihution  for  the  area  ovct  w  Inch  the 
pulse  e.vpands  remains  (iaus.-ian  throughout,  the  shajie  of 
the  signal  -S'/.i  resulting  from  all  scatterer.s  i-  given  hv  the 
convolution  of  function.-  describing  tlie  transmilli.'d  jiul.-c 
shape,  the  antenna  gam  pattern,  and  the  sui  lace  height 
distribution  I  Dnnkwaler,  1991 ).  Thi.s  .signal  i-  the  pulse 
echo  waveform  •  Figure  7-2 1.  The  half-|K)wer  point  is  at 
delay  time  /o,  which  give.-  (he  range  to  mean  .sea  Itwti. 
Maximum  power  occurs  at  the  delay  time  at  which  the 
maximum  footprint  tirea  i.s  achieved.  At  later  delay  time.-, 
power  drops  off  due  to  the  antenna  gam  pattern.  Over  the 
ocean,  the  slope  of  the  leading  edge  tor  ramp)  oi  the  w  ave- 
form  is  inversely  proportional  to  SWH.  and  the  peak  power 
i.s  related  to  wind  speed.  The  characteristic  shape  of  lhi.s 
signal  only  holds  for  reflecting  .surfaces  of  meter-scale 
roughnes,-  for  which  specular  point  .scattering  can  be  as¬ 
sumed. 

The  digitally  sampled  "■a'.'cform  has  a  ’•■due  -at  60  'for 
Seasaticonliguous  time  intervals,  centered,  ifthe  tracker  is 
operating  correctly,  on  the  half-power  point.  The  time 
intervals,  or  range  gates,  each  correspond  to  a  two-way 
travel  time  resolution  of  3.125  nano.seconds,  which  is 
equivalent  to  a  range  resolution  of  about  1  motor,  i  Better 
precision  in  the  range  to  nadir  or  height  measurement  is 
obtained  through  averaging).  The  powerintherangegates 


Time 


Fig.  7-2.  The  pulse  echo  for  pulse-limiteH  .sU'rneiry  and  a  Gaui,-inn 
surface  heiglil  di.stribution.  The  half-power  point  at  /„  correspond? 
to  the  delay  of  return.?  from  the  mean  sc, a  level 


is  the  result  of  backscatter  from  range  rings  within  the 
altimeter  footprint  (Figure  7-3).  The  area  covered  by  each 
range  ring  is  equal  to  that  of  the  PLF 

/Ipif-.  Tirpif  =:  )xc//r  (1) 

where  Tpif  is  the  radius  of  the  PLF,  His  satellite  height  above 
the  surface,  and  ris  compressed  pulse  length  (neglecting, 
for  simplicity,  the  Earth’s  cuiwature).  With  this  expression 
for  area,  it  is  possible  to  use  the  radar  equation  to  arrive  at 
an  expression  for  the  power  received  from  any  range  ring 


P 


range  ring 


P, - ;; — — o 

4(47T)2//3 


(2.) 


where  Latm  is  atmospheric  attenuation  loss.  Using  a  form 
of  this  equation,  backscatter  can  be  retrieved  from  wave¬ 
form  range  gate  samples.  For  Geosat,  the  formula  is 

=  lOlogjQ  Pi  +  AGO  +  K  +  SOlogjQ  Lattn  (3) 

where  Pjis  the  power  in  range  ring  i  expressed  incounts,  Hq 
is  a  reference  height,  AGC  is  the  attenuator  setting  (de- 


Fig.  7-3.  Altimeter  footprint  range  rings  (not  to  .scalei.  Each  ring 
corresponds  to  a  range  gate  in  the  digitally  sampled  waveform. 
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scribed  in  Section  7.2.4).  and  K  i.s  a  con-stant  dependent  on 
system  parameters  lUlander.  1  987u |.  'I'he  peak  back.i^cai- 
ter  usually  appears  in  one  ol’the  first  few  l  unge  ring.--  for  .sea 
ice  and  is  generally  a.s.sumed  to  be  the  back.scatter  from 
nadir.  When  altimeter  back.scatler  rneasuretnent.s  are 
referred  to  in  what  follows,  it  is  u.suallv  the  peak  back.scat- 
ter  that  is  meant. 

While  the  PLF  corresponds  to  that  part  ofthe  waveform 
over  which  the  maximum  power  is  reached,  the  range 
window-limited  footprint  corre.sponds  to  the  area  nuqiped 
out  by  the  pulse  at  the  end  of  the  time  corresponding  to  the 
last  range  gate.  For  Seasat,  the  diameter  of  the  range 
window-limited  footprint  was  about  1 0  km,  and  the  an.gie  of 
incidence  at  the  last  range  ring  was  about  0.3  a 

7.2.3  Wauefonn  Areraging 

In  the  case  of  noncoherent  scattering,  the  phase  relat  ion  - 
ships  of  returns  from  all  facets  or  specular  points  arc- 
random.  The  vector  sum  of  amplitude  and  pha.se  for 
reflections  within  a  range  ring  is  different  for  each  range- 
ring,  which  lends  noise  or  fading  characteristics  to  indi¬ 
vidual  pulse  echoes.  Therefore,  waveform.s  are  averaged  to 
reduce  fluctuations  and  to  provide  a  good  estimate  of  heigh  t 
(given  by  the  half-power  point)  and  backscatter,  Lacn 
waveform  must  be  a  statistically  independent  .sample.  The 
pulse  repetition  frequency  ( PRF  i  ncces.sarv-  to  avoid  coitc- 
lation  between  pulse  echoes  depend.s  on  the  surface 
roughness  and  orbital  velocity,  Sea.sat  used  a  PRF  of 
1020  Hz.  Fifty  echoes  were  averaged,  during  which  tirne 
the  satellite  nadir  track  moved  about  330  m.  Foraltirrudrv 
over  the  ocean,  averaging  50  echoes  is  a  good  compromise 
that  provides  enough  .samples  for  an  accurate  power  esti¬ 
mate  while  not  degrading  the  along-track, spatial  resohil  ion 
of  the  measurement  beyond  acceptable  bounds,  given  the 
typical  scale  of  changes  in  dynamic  sea  .surface  height,  winu 
speed,  and  SWH.  When  surface  charactcristic.s  an.  changing 
over  the  integration  time  (us  is  often  the  ca.'C  ovw  in-  or 
when  pulse  echoes  are  dominated  by  one  or  .--everal  .-trong 
scatterers,  waveforms  that  result  from  averaging  will  not 
truly  represent  the  surface. 

7.2.  i  The  Adaptive  Tracker  and  Prov’.^sing  Lonjjs 

Accurate  estimates  of  height  and  other  parameter.'  dc 
rived  from  waveforms  depend  on  accurate  superposition  of 
pulse  echoes  for  averaging.  This,  in  turn,  depends  on  proper 
timingofthe  deramping  pulse.  The  following  description  ■! 
the  method  by  which  this  is  achieved  is  liased  on  the  at 

adapLv  ■  tracking  unit,  but  Geosat  and  KRiS-I  trackers  ;u c 
similar  in  concept.  The  adaptive  trai  l,,,  i  i  i  avh.  •o'-ttic 

through  an  automatic  gain  control  (.\Gt ' '  loop  an<i  trai  kine 
(or  range  and  range  rate  i  lonp.s.  Tin-  AGG  conqionent  of  t  he 
tracking  system  constantly  adjust, s  the  ,  v 

Loacliieve  It  le  maxim  inn  sutnai-io-noise  ratio.  IVnvi  r  wiildn 
the  (averaged.!  waveform  is  computed  accordin:;  to 
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AGC\..,u-  =  '4 1 

i-i 

where  i  is  range  gate  number,  P  is  the  power  iexpress.  il 
as  a  count  value)  in  each  gate,  and  Sc  is  a  scale  I'actor,  Fig¬ 
ure  7-4ia).  An  attenuation  value  is  computed  according  to 
the  recursion  ''elation 

AGC  If!  -r  1 )  =  AGC  in)  +  aAdtJC  < 5 1 

where  AAGC  is  the  dilTerence  between  and  some 

reference,  and  a  is  some  constant  1 0.02hri  for  Seasat  >.  The 
AGC  loop  works  to  keep  the  average  po'»ver  in  t  he  waveform 
constant.  It  is  updated  at  a  rate  of ’20  Hz  i  every  50  pulses ). 
and  the  new  value  is  applied  to  the  next  50  piil.se  echoes. 
Because  of  the  factor  a  in  the  recursion  relat  ion.  an  adju.st- 
ment  to  the  value  of  AGC  necessitated  by  a  change  in 
backscatter  from  the  surface  takes  place  over  O.S  secimds, 
during  which  time  the  satellite  moves  5.3  km.  While  thi- 
parameter  AGC  is  a  good  estimator  of  backscatter  over  the 
ocean,  this  is  not  the  case  over  ic(>,  where  peak  power  can 
change  an  order  of  magnitude  in  a  much  .shorter  di.st  'ince. 
Therefore,  to  retrieve  high-frequency  \';ir(iibi!ily  in  back¬ 
scatter  it  is  neces.sary  to  take  into  account  the  count  vaiue.s 
in  each  waveform,  along  with  the  AGC  value  that  was  the 
attenuator  setting  during  the  waveform  integration  time. 
In  such  cases  where  backscatter  is  changing  quickly,  A(»C 

(a) 


Fig.  7--1.  Tile  altimeterrange  windov.  and  gates.  (alOvertheocean, 
I  he  !.>•■■  -  I  111  -hr- a  a  ...q  rdi'i  he  power  in  *he  waveform. 

The  half-power  point  occurs  at  time  f,,  if  the  tracker  is  operating 
correctly.  '  hi  Early  ( ft  >,  middle  <  M  i,  and  late '  h'  gates  in  the  center 
of  the  waveform  detennine  the  slop-  of  the  ramp  ( proport inn,al  to 
SWH I.  The  half-power  point  .should  be  in  gate  M. 


Is  overeslimaliiig  or  uiidei'i  si  mi.iti!:;'  power,  .uo!  !f-  p: 
nal-to-iiuise  I'.ilm  v.  il!  not  b,.  ..jniniol 

In  Figure  7-1  itypic.il  oi  oio.ui  -  n'liiil.-  .  li!.  isuiii-  t,: 
.4GC.,,j,.  is  about  t-iju.ij  III  hah  tin  p,,n,-r  ju  lie  w,;\  i  ion), 
The  power  in  t  he  ct-nti-r  of  I  hi-  r.in.e  .Miidow  i-  oa-.i-u!,,! 
Ill  a  middle  gate  .51,  wliictl  Miriietiisie'  >.  . tiled  gal,  S'  a 
Figure  7  -l  b-  Tile  ll-ack,-!  .seek-  to  triinionee  Ig,  I 
error  \r,  where 

V-nl(,(' 

by  sliilting  tlu'  range  uindov.  in  time  to  Keep  tia  halt 
power  point  in  .\1.  I’he  w  idtlud  gate  M  depend;  i-n  1 1  ,.  •  i<  .pt 
ofebe  leadmgedgei.d'th.e  wa',  etonn.  .iiid  tiiei  ■  ti.ie  on  .sW  M 
The  tracker  adapt.--  the  w  idt  li  id  M  m  range  gate-  !>a  ed  oi: 
the  dilVerence  between  pniwer  in  an  earl\  ,;,iti  h  -i'  the 
bottom  oft  hi  ramp  ;ind  a  later  ;:.iU  1.  .it  '  he  tope!  tla  -  t.ngp 
Fora  -significant  wave  hei-.dit  ofb  th-m  2  m,  ::.tie  M  g 
.single  range  g-ite  wide  i  t 'helton  et  a!  ,  I'.tstl 

fhe  range  error  .\r  i-  ii-ed  to  pi.-iIh  t  tia  !  to  tie 

surtace'lind  therefore  thie  delay  tune  .d  the  e,  sp,  i  o'  tia 
range  vmdo'w  ■  loi  .-uiisi  qm-nt  pui.-e-  .m  ordnig  to 
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when-  .\f  i.s  the  range  error  lor  juil:  e  c  .  a  g-.  the  i r.u'k'  t  loni 
constant,  M  is  tlie  update  interval  of  the  I'Kip,  and 
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i.s  the  range  rale.  Time  con.stants  a  tind  b  c.ui  be  adpi-ted 
to  incrr.'ise  the  ;tgility  of  fhe  tr;icker  -  at  the  price  ot  in¬ 
creased  noi.sf'.  Din'erent  design.s  for  trticking  ocean  .inii 
nonocean  surface--  have  been  explored  IMclnfcre  et  ak. 
1983;  Kapley  et  al..  19S-3;  I’lander.  !‘,)87t:!.  The  FH.''  ] 
;iUimeter  has  both  an  oci'an  mode,  which  use,-  a  .Suhopt  ima! 
Maximum  Likelihood  K,-tim:itor  iSMl.K:  tracker,  and  an 
ice  mode,  which  has  a  spht-gtili'  tr.-icker  unlike  the  Se.i-at 
Geo.sat  tracker  I  Gnffith.set  ak.  1  iiskl.  The  i<  e  mode  tracker 
has  a  wider  range  window  for  better  track  in;:  o\  I  r  uruiulatmg 
continental  ice  shelves.  For  .“-ea  let',  howevi-r.  it  g-  not  h<  ;ght 
vanation.s  confounding  the  agihty  oft  he  tracker,  luit  i  other 
backscatter  variations  within  the  hiotprint.  whicii  intro 
duce  errors  oxacerhateti  by  interaction  between  the  j're 
cessing  loops. 

7.2.5  firror.s  Common  (.Irrr  h  r 

Figure 7-.5  illustrates  how  the  rapid  n-e  and  f;dl  of.soa  ice 
waveforms  eonioun'i.)  llie  irackei  .--eiei.  iioo  of  go. L,  7vl , 
and  L.  Because  of  the  steepness  of  the  ramp  and  the  rafiid 
fall-off  from  the  peak  value  in  the  return  waveform  in 
Figure  r-.aia  i,  the  power  in  M  is  not  half  the  avertige  power 
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in  the  waveform.  Power  within  the  AGC  gate  is  much  less 
than  that  within  the  middle  gate.  In  extreme  cases,  this 
results  in  a  large  negative  height  error,  causing  the  tracker 
to  shift  the  range  window  .so  that  .subsequent  pulses  appear 
later  in  the  range  window.  Figure  7-5(b.).  On  the  next 
update  cycle,  the  middle  gate  measures  the  noise  before  the 
ramp,  resulting  in  a  large  positive  height  error.  The  tracker 
corrects  for  this  until  the  peak  is  in  the  middle  gate  again 
when  the  error  repeats.  This  error,  called  tracker  oscilla¬ 
tion,  is  evident  in  data  that  exhibits  height  excursions. 
Tracker  oscillation  can  sometimes  be  corrected  by  relracking 
the  waveforms.  Ketracking  involves  (1 )  appl>ing  an  algo¬ 
rithm  to  the  waveform  record  that  finds  the  point  on  the 
waveform  that  should  corresponds  to  return  from  the  sur¬ 
face  at  nadir,  and  (2 )  using  the  position  of  that  point  relative 
to  the  center  of  the  range  window  to  correct  the  height 
estimate  based  on  the  range  window. 

The  range  and  range  rate  loops,  like  the  AGC  loop,  w’ere 
updated  at  20  Hz  (every  50  pulses)  on  Geosat  and  Seasat. 
Pulse  echoes  acquired  with  the  same  loop  settings  were 
averaged  to  make  waveforms.  However,  waveforms  were 
telemetered  at  10  Hz,  so  the  telemetered  waveform  is  the 
average  of  two  50-pulse  mean  waveforms  from  two 
tracking  loop  updates.  If  the  position  of  the  range 
window  is  changing  rapidly  (due  to  changing  topography  or 
tracker  oscillation),  the  resulting  telemetered  waveform 
will  be  a  blurred  or  sometimes  double-peaked  composite, 
Figure  7-6(a).  Errors  resulting  from  the  incorrect 
superposition  of  waveforms  for  averaging  are  termed  te¬ 
lemetry  summing  errors  [Rapley  et  al.,  19871,  ERS-1 
telemeters  waveforms  at  the  tracker  loop  update  rate)  20  Hz), 
thereby  avoiding  the  problem  that  appeared  in  Seasat  and 
Geosat  data. 

(a) 


Range  Window 

1 

(b) 


Fig.  7-5.  (a)  The  rarnpofa  typically  peaked  waveform  from  sea  ice 
can  appear  entirely  in  gate  M.  (b)  This  causes  the  tracker  to  shift 
the  range  window  earlier  in  time  for  the  next  update  cycle. 


Tracker  o.scillation  caiuiut  be  coinpletelv  ehrniiialed  by 
retracking  if  telemetry  summing  error.s  are  pre.seiil,  Fig¬ 
ure  7-f>(bi.  However,  the  probability  that  .such  error.s  exi.sl 
can  be  estimated  L>  examiiiing  llie  range  rate  parameter, 
which  will  peak  when  the  range  window  i.-  clianging  at  the 
greatest  rale  and  the  probability  of  blurring  i.s  highest 
[Laxon,  1989!.  The  combined  effect  of  tracker  o.sf  illalion 
and  telemetry  summing,  w  hen  it  produces  a  .sudden  .shill  in 
range  and  blurred  or  double-peaked  waveforni.s,  is  some 
time.s  referred  to  a.s  height  glitch, 

VVhen  the  AGC  loop  set.s  the  .AGC'  value  Kki  low  to  ac¬ 
commodate  power  in  the  peaked  port  ion  of  an  ice  waveform, 
saturation  occurs.  Saturation  was  common  in  .Sea.sat  data 
iWingham  and  Rapley,  I987|.  It  is  often  a.'.,si>cialed  with 
frequency  splatter.  Figure  7-8(0.  where  power  from  the 
saturated  gate  appears  in  earlier  and  later  gate  -  Vnother 
artifact  often  present  in  ice  data  i.s  .snagging  or  "off  ra  ngi  ng ,  ■ 
This  occurs  when  a  bright  feature  al  nadir  is  follow  ed  by  ( he 
tracker  as  it  recedes,  until  it  is  sufficiently  far  from  nadir  for 
the  true  nadir  return  to  once  again  dominate.  Often  ih..s 
occurs  at  a  transition  from  pack  to  fast  ice.  or  from  thin  ice 
to  multiyear  (MY)  ice  (Laxon,  1989). 

Antenna  off-pointing  was  a  serious  problem  with  Geo.sat , 
Occasionally,  the  off-pointing  was  severe  enough  to  cause 
loss  of  lock,  a  situation  where  the  waveform  disappears 
from  the  range  window  entirely  for  a  few  telemetry  fraine.s. 
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Fig.  7-6.  Common  errors  in  altimetry  over  sea  ice.  (a)  An  example 
of  telemetry  summing  cau.sing  double  peaking  in  a  Seasat  wave¬ 
form  over  ice  I  Ulander,  1987b  |.  (b)  Sea.sat  data  .showing oscillation 
in  the  onboard  height  measurement.  Retracking  does  not  remove 
the  oscillation  due  to  telemetry  summing  errors  I  Laxon,  19891. 
(c)  A  saturated  Seasat  waveform  [Ulander,  1987b|, 
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Over  the  ocean,  off-pointing  can  be  identified  by  the  atten¬ 
dant  rise  in  power  in  later  gates,  and  a  correction  can  be 
applied.  Over  sea  ice,  unpredictable  wavefonn  shapes 
make  this  impossible,  resulting  in  altered  waveform  shapes 
and  unidentifiable  errors  in  height  and  buckscatter  mea¬ 
surements.  The  stabilizer  used  on  ERS-l  should  prevent 
off-pointing  from  being  a  serious  problem. 

Once  the  altimeter  record  over  ice  has  been  retracked 
and  edited,  the  backscatter  coefficient  for  any  range  gate 
can  be  found  from  the  wavetorm  using  Equation  (3)  with  the 
telemetered  w'aveform  range  gate  count  values  and  the 
value  of  the  attenuator  setting  (AGO).  However,  it  is 
important  to  remember  that  the  altimeter  moves  aboui 
0.66  km  between  telemetered  waveforms  (for  Geosat  and 
Seasat).  If  the  statistics  of  the  surface  are  not  stationary 
over  that  distance,  the  backscatter  estimates  may  be  mis¬ 
leading. 

7.3  Sea  Ice  ScArrERLNG  Models  at  Normal  I.nctde.n'ce 
7.3. 1  Introduction  to  Altimeter  Backscatter  Models 

The  distinguishing  characteristics  of  ice  waveforms  are 
that  they  are  often  narrow-peaked  in  shape,  taper  off 
quickly  with  incidence  angle,  have  high  backscatter  coeffi¬ 
cients  at  nadir,  and  are  widely  variable.  Average  annual 
values  of  Seasat  backscatter  over  Arctic  and  Antarctic  sea 
ice  ranged  from  approximately  0  to  45  dB,  with  a  mean 
greater  than  20  dB  [Laxon,  1989j.  In  contrast,  backscatter 
over  the  ocean  has  a  range  of  between  8  and  15  dB, 
depending  on  sea  state  [Chelton  and  McCabe,  1 985].  On  the 
basis  of  dielectric  characteristics  alone,  one  would  expect 
ocean  backscatter  to  be  greater  than  that  for  ice,  and  indeed 
for  surface  measurements  of  backscatter  at  nadir  this  is 
usually  the  case.  The  dynamic  range  of  altimeter  measure¬ 
ments  over  ice  cannot  be  explained  by  variability  in  the 
dielectric  characteristics  of  ice  alone.  Since  backscatter  for 
airborne  and  satellite  altimeters  is  dominated  by  surface 
rather  than  volume  scattering,  the  nature  of  altimeter 
pulse  echoes  must  be  primarily  determined  by  surface 
roughness.  Specular  or  coherent  reflection  from  smooth 
surfaces  near  nadir  may  produce  the  peaked  component  of 
waveforms,  while  diffuse  reflection  from  rough  surfaces 
may  account  for  power  in  later,  off-nadir  gates.  But  while 
sea  ice  or  ocean  surfaces  can  in  certain  circumstances  be 
extremely  smooth,  the  question  remains  as  to  whether 
these  surfaces  are  smooth  over  large  enough  areas  to 
contribute  coherent  returns.  Models  that  invoke  only 
noncoherent  (diffuse)  reflection,  on  the  other  hand,  may 
have  difficulty  describing  the  high  backscatter  from  nadir 
and  narrow  peaked  shape  often  .oeen  over  ice. 

The  physical  characteristics  of  sea  ice  are  quite  variable, 
which  compounds  the  modeling  problem.  Sea  ice  consists  of 
both  level  and  deformed  ice  of  different  ages,  often  with 
superimposed  dry  or  wet  snow  layers.  Water  openings  of 
different  sizes  and  with  variable  wind  conditions  are  often 
present.  The  different  layers  and  interfacesofthe  snow,  ice, 
and  water  may  all  contribute  to  the  .scattering  of  electro¬ 


magnetic  waves.  Depending  on  the  frequency  u.sed,  how¬ 
ever,  It  i.s  often  po.ssible  to  neglect  particular  scaltenng 
terms  and  thus  simplify  the  .scattering  computation  For 
example,  the  surface  will  appear  rough  and  result  m  the 
donunance  of  diffuse  scattering  if  the  .surface  height  v  ana- 
tions  are  large  compared  with  the  radar  wavelength.  The 
dependence  of  surface  roughnes.s  upon  radar  wavelength  i.s 
expressed  in  parameters  fen  and  kl,  where  A*  =  2tja.  a  i.s  the 
standard  deviation  of  .surface  height,  and  I  i.s  the  correlation 
length  of  surface  roughnes.s.  The.se  pararneter.s  determine 
the  range  of  validity  for  coi,.monly  used  model.s.  In  the 
following,  we  consider  surface  scattering  clo.se  to  normal 
incidence,  The  emphasis  i.s  on  radar  wavelength.s  i  n  the  Ku- 
ba  .d  <1 .7  to  2.4  cmi.  since  most  radar  altimeters  u.sed  for 
geophysical  remote  .sensing  operate  in  this  band. 

7.3.2  Backscatter  Front  a  Stationary  (iuus.stan  Random 
Surface 

This  section  concerns  backscatter  from  level  and  snow  ■ 
free  sea  ice.  The  ice  is  modeled  as  a  homogeneous  medium 
wdth  constant  dielectric  constant.  Its  surface  has  small- 
scale  height  variations  represented  by  a  stationary  Gau.ssian 
process.  The  electromagnetic  scattering  from  .such  a  sur¬ 
face  has  been  studied  extensively  in  the  past .  but  a  general 
closed-form  solution  has  not  yet  been  found.  A  recently 
developed  model  for  a  perfectly  conducting  [Fung  and  Pan. 
1987)  and  dielectric  [Chen  and  Fung,  1991]  surface  ha.s 
promise,  bat  its  range  of  validity  has  not  yet  been  deter¬ 
mined.  High-  and  low-frequency  solutions,  however,  arc 
readily  available  and  are  sufficient  for  the  pre.sent  discus¬ 
sion. 

7.3.2. 1.  Specular  backscatter  from  a  perfectly  smooth 
surface.  A  surface  is  usually  considered  smooth  if  it  meets 
the  Rayleigh  criterion.  Fo'  Ku-band.  a  must  be  less  than 
abo  ut  0.3  cm.  The  scattering  pattern  for  such  a  surface  will 
show  a  large  phase-coherent  component  of  backscatter  in 
the  specular  (nadir)  direction  and  a  much  smaller  diffuse 
component  at  other  incidence  angles.  If  a  surface  is  per¬ 
fectly  smooth  and  planar,  scattering  is  confined  to  the 
specular  direction,  with  no  diffuse  component,  and  is 
equivalent  to  ray  reflection.  Such  a  surface  will  produce 
true  specular  reflection  if  it  is  perfectly  smooth  over  at  least 
two  Fresnel  zones.  The  edge  of  each  circular  Fresnel  zone 
is  defined  where  the  phase  of  the  reflected  wave  has  changed 
by  0.571  compared  with  a  reflection  from  the  nadir  point. 
The  diameter  of  the  first  zone  is  given  by 

=  (9) 

or  about  180  m  for  an  ERS-l  class  altimeter.  Although 
possible  in  theory,  it  is  unlikely  that  the  rms  height  of  the 
surface  is  perfectly  smooth,  or  even  smooth  to  within  a  few 
millimeters,  over  an  area  of  several  hundred  square  meters. 

It  is  not  convenient  to  describe  true  specular  reflection  in 
terms  of  a  backscatter  coefficient.  However,  it  is  possible  to 
derive  an  effective  backscatter  coefficient  that  is  dependent 
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on  the  radar  system.  For  example,  the  backscatter  coeffi¬ 
cient  of  an  infinitely  large  dielectric  mirror  observed  by  a 
pulse-limited  altimetei  is  given  by  [Ulander,  1987al 

cr’eoh  =  l/?(0)|2—  (10; 

cr 

where  i?(0)  is  the  Fresnel  reflection  coefficient  Rid)  evalu¬ 
ated  at  0°.  The  subscript  indicates  that  reflection  is.  of 
course,  coherent.  The  theoretical  maximum  return  for  an 
ERS-1  class  altimeter,  then,  is  about  48  dB,  assuming  an  ice 
reflection  coefficient  of  about  -1 1  dB. 

7. 3. 2. 2.  Coherent  backscatter  from  a  slightly  rough 
planar  surface.  When  a  small  rms  surface  roughness  is 
superimposed  on  a  flat,  planar  reflecting  surface,  an  ex¬ 
pression  may  be  derived  for  the  coherent  backscattering 
coefficient,  acoh*  which  accounts  for  the  sphericity  of  the 
wavefront  and  the  radar  system  impulse  response  [Fung 
and  Eom,  1983]: 


iRieV 

tTcoh  = - ; - exp 

Ppir 
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(11) 


Equation  (11)  assumes  a  Gaussian  beam  describud  by  the 
one-sided  pulse-limited  beam  width  where 
At  normal  incidence  this  simplifies  to  [Ulander  and 
Carlstrom,  1991 J: 


crVoh  =  I  ^<0)  1 2  —  exp  d  2) 

CT 

Figure  7-7  illustrates  the  coherent  backscatter  function 
derived  using  Equation  (11)  with  suitable  parameters  for 
the  ERS-1  altimeter  and  the  Rutherford  Appleton  Labora¬ 
tory  (RAL)  airborne  altimeter  [Drinkwater,  1987]  respec¬ 
tively.  For  the  satellite  case,  the  fall-offin  <7°coh  extremely 
rapid  for  all  surface  roughnesses.  The  maximum  backscat¬ 
ter  values  of  between  40  and  48  dB  occur  for  slightly  rough 
surfaces  typifying  new  saline  ice  (i.e.,  f‘*'=4.0  +  7O.4).  In 
contrast,  the  RAL  altimeter,  with  a  pulse-limited  beam 
width  ofO.GS",  gives  a  much  wider  response  function  that  is 
lower  in  peak  power  at  all  roughnesses.  These  results 
demonstrate  the  role  that  the  compressed  pulse  length  and 
altitude  of  the  altimeter  play  in  dictating  the  coherent 
response  function  of  a  slightly  rough  surface.  In  short,  the 
system  impulse  response  is  the  most  important  factor  in 
determining  the  decay  of  the  coherent  scattering  coefficient 
with  incidence  angle.  Furthermore,  the  coherent  power  is 
sharply  reduced  as  the  surface  roughness  increases  only 
slightly. 

7. 3.2. 3.  Noncoherent  backscatter  from  a  rough  surface 
where  kl  is  large,  but  ko  is  small.  Field  measurements  of  sea 
ice  surface  roughness  typically  show  that  the  rms  height  a 
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Fig.  7-7.  Coherent  backscjitter  versus  incidence  angle  for  .spiice  and 
airborne  platforms. 


is  0.1  to  1.0  cm  and  the  height  correlation  function  is 
exponential  or  Gaussian  with  a  correlation  length  I  of  2  to 
10  cm  [Ulander,  1991).  For  radar  wavxdengths  in  the  Ku- 
band,  this  ice  surface  is  in  the  intermediate  roughness 
range.  Note  that  the  correlation  length  is  generally  larger 
than  the  radar  wavelength.  The  KirchholT approximation 
can  therefore  be  used  to  evaluate  the  scattering  integrals 
(Beckmann  and  Spizzichino.  1963).  This  approach  uses  the 
diffraction  integral  for  the  scattered  field  and  the  tangent- 
plane  approximation  to  compute  the  surface  currents.  The 
latter  requires  that  the  average  radius  of  curvature  of 
scattering  elements  be  large  compared  to  the  wavelength. 
The  rms  surface  slope  must  also  be  less  than  0.25  radian.s 
in  order  to  allow  computation  of  an  approximate  (scalar.i 
solution  for  dielectric  surfaces.  With  these  restrictions,  the 
noncoherent  backscatter  coefficient  cr'^ori  given  by  [  Ulaby 
etal.,1982] 


c/*{ 0)  =  Ir(  0)1^  2k  ^cos'^0  exp(-g)  ^  W "  ( 2k  sin0,  o) 


(13) 


whereg  =  t2kacos6)^  and  the  first-order  slope  term  that  is 
important  for  larger  incidence  angles  has  been  excl  uded .  I  n 
Equation  (13),  the  roughness  spectrum  Ky)  is  re¬ 

lated  to  the  power  of  the  correlation  function  defined  by 
(Fung  and  Pan,  1 987 [ 


W"  (A„ K,)  =  -U"  f”  p--  ( I,  f )  exp  i-jKy^  -jKyC)  dldf  ( 1 4 ) 

J  - 

where  pt^.C)  is  the  two-dimensional  normalized  height 
autocorrelation  function.  Analytical  expressions  for  the 
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iDte^Tal  iji  Eqiuuion  i !  4),  when  the  surface  is  isulropic  and 
has  a  Gaussian  or  exponential  orrolation  function,  can  be 
found  in  Ulaby  et  al.  [19821  and  Kim  ot  al.  |1985|,  re.spec- 
tively. 

7. 3. 2.4.  Noncoherent  backseat  ter  from  a  rough  surface 
whc’-e  kl  and  ka  arc  large.  In  the  Kirchhoff  or  phy.sical 
optics  theory,  when  rms  height  is  much  larger  than  the 
radar  wavelength  it  is  possible  to  evaluate  the  diffractioti 
integral  usinga  stationary  phase  approximation  [Beckmann 
and  Spizzichino,  19631.  This  coiTO.sponds  to  the  geometric 
optics  limit,  in  which  the  .surface  is  assumed  to  be  made  up 
of  tangent  planes  or  specular  facets  that  are  inf  initely  large 
with  respect  to  a  wavelength.  Thus  it  is  often  known  as 
specular  point  theory.  For  such  rough  surfaces,  backscatter 
is  completely  noncoherent  and  the  backscatter  coeftlcient  is 
given  by  [Ulaby  et  al..  1982] 

(rt  0  )  =  exp  (-  feb-ifi j  ,15, 

.s-cos-0  '  ' 

where  s-  is  rms  slope  in  radians. 

The  scalar  and  tangent-plane  approximations.  Equa¬ 
tions  (13)  and  (15),  respectively,  both  give  the  mean 
noncoherent  backscatter  from  a  stationary  Gaussian  ran¬ 
dom  surface.  Since  backscatter  is  being  summed 
noncoherently,  if  an  altimeter  were  to  .scan  such  a  rough 
surface  the  returned  power  would  exhibit  statistical  fluc¬ 
tuations  or  signal  fadingidescribed  in  Section  7.2.3).  The 
return  power  is  then  distributed  according  to  a  negative- 
exponential  or  Rayleigh  distribution,  provided  that  the 
following  conditions  are  satisfied; 

(1)  Many  individual  scattering  elements  contribute  to  the 
returned  power,  and  no  single  one  dominates. 

(2)  The  amplitude  and  phase  of  returns  from  the  scatter- 
ers  are  independent,  and  the  phase  distribution  is 
uniform. 

7.3.3  Backscatter  From  a  Heterogeneous  Random  Surface 

Equations  ( 1 1 )  and  ( 1 2)  show  that  the  coherent  return  is 
a  highly  nonlinear  function  of  the  rms  surface  height.  When 
the  rms  height  changes  only  from  0.2  to  0.4  cm,  the  coherent 
component  is  reduced  by  17  dB.  This  suggests  that  the 
normal-incidence  response  is  dominated  by  the  smoother 
parts  of  the  surface  for  which  the  rms  height  is  a  tenth  of  a 
wavelength  or  less.  The  first  model  to  recognize  this  fact 
was  that  proposed  by  Brown  [19821.  The  model  assumes 
that  the  .scattering  originates  from  an  ensemble  of  circular, 
horizontal,  and  flat  dielectric  patches  near  nadir.  The 
patches  area.s.sumed  to  cover  a  .small  fraction  ofthesurface, 
to  be  located  at  the  same  height,  and  to  be  randomly 
distributed  in  the  horizontal  plane.  Brown  shows  that  both 
a  coherent  and  a  noncoherent  component  result  from  this 
model.  The  former  is  caused  by  the  patches  being  located  at 
the  same  vertical  level,  to  within  a  few  mm,  and  is  propor¬ 


tional  to  the  numbi'i'  of  patche,--  .-quarod.  'J'lu  !atu  s'  i.- 
caused  by  the  random  horizontal  d is!  rib' a  lonofslu-p.itt  li*-.- 
and  i.s  directly  pioporiional  icj  the  number  of  pairlie.-. 
.Although  the  model  ri  co,.,.  z(a:  tin-  im])or lance  <4  !  in- 
tia!  vanabiiity  of  the  surface  chai-actiaa.--ta.-,'..  il  nihei'eu!  i;. 
a.ssume.s  that  rellection.--  fi-oin  the  palclues  add  colu  ia  tuh  , 
Tin  .sas.sumption  is  probably  unre.ilislie  in  [iratlice,  at  lea  -t 
for  .spaceborne  .sy.stems,  due  to  the  large  llr.st  Fre.'iiel  zone 
Another  model  ha.s  been  propo.sed  which  does  n<  it  n-qu  i  re 
the  fiat  patches  to  add  coherent  ly  [  Ulande!'  and  (..'arlsl  rum . 
1991 1.  The  back.scattered  powei'  i.s  in.st<-ad  a.s.sunu  d  to  be 
dominated  by  refiections  from  horizontal  and  ilat  patehe.'- 
that  add  noncoherently.  narrow-peaked  waveform  nia\ 
occur  ifthe  palche.s  are  .sufficiently  large.  This  a.ssumpUon 
i.s  more  realistic,  but  field  rnea.surement.s  are  presently 
lacking  to  verify  it.  Using  phy.sical  optics,  the  radar  cro.-.- 
section  from  a  single  circular  and  fiat  dielectric  patch  is 
given  by 

aid}  .4  R  ( ( 1 )  I  ^  ^  j  IX  p|  I  '  J  (i  i 

where  A  i.s  the  patch  area.  I)  is  the  patch  diameter,  and  t  he 
angular  respon.se  has  been  approximated  In  a  Gau.-i>ian 
function  with  ---  U..3  hi).  For  a  collection  of  .sparsely 
distributed  and  horizontal  patches,  we  then  obtain  the 
backscatteringcoefficient  according  to 

where  F  i.s  the  fractional  coverage  of  flat  patches. 

Equation  (ITi  does  not  take  the  altimeter  sy.sli. m  im¬ 
pulse  response  into  accouni.  Wht*n  the  patch  di.imeter  is 
sufficiently  large  to  produce  a  narrow-peaked  wavcl'orm, 
the  effective  baclvscattering  coefficient  is  given  by 

a"  =  1-  i’  'dr  =  d.AFR  0'^^  a8i 

r„  Uli  cr 

Typical  values  for  the  narrow-peaked  backsciitter  coeffi¬ 
cient,  as  measured  by  spaceborne  altimeters,  are  20  to 
40  dB,  corresponding  to  FA).  2(1  to  16'.;  according  to  Equa¬ 
tion  (18.1.  In  this  class  of  .surface,  the  coherent  component 
can  exhibit  signal  fading  as  the  number  of  palche.s  or 
reflecting  facets  within  the  first  Fresnel  zone  vanes  from 
pulse  to  pulse.  Although  the  terms  specular  and  coherent 
reflection  are  often  used  interchangeably,  true  specular 
reflection  originates  from  a  single  infinitely  large  (when 
compared  with  the  FLfT  mirror- 1  ike  facet  and  doc.s  not  show 
fluctuations  in  amplitude  iind  pha.se  from  pulse  to  pulse 
U  is  also  pos.sibic  to  modify  the  Kirchhoff  models  accord¬ 
ing  to  Equations  (13)  :ind  <1  5)  to  inc  lude  the  effects  of  a 
heterogeneous  surface  and  to  feoduce  a  nairow-peaked 
waveform.  We  here  give  the  result  for  the  geometric-optic.s 
model 


wWch  only  differs  by  a  factor  0.9  from  Equation  (18).  Ttds 
shows  that  the  geometrical  optics  and  flat  patch  model 
predictions  are  indeed  very  similar,  although  the  scattering 
mechanisms  assumed  are  very  different.  The  approach 
taken  by  Ulander  and  Carlstrbm  [1991 1  is  more  realistic  in 
that  it  accounts  for  height  distributions  or  surface  tilting  of 
patches  and  is  not  restricted  by  the  approximation  of 
scattering  by  individual  facets  to  the  scattering  from  smooth 
specular  points.  But  while  the  assumptions  of  the  other 
models  are  simplistic,  the  notion  that  scattering  from  many 
individual  surface  element.s  conforms  reasonably  closely  to 
that  from  a  specular  surface  appears  not  too  far  from  ihe 
truth. 

7.3.4  Discussion 

In  order  to  contrast  normal  incidence  coherent  and 
noncoherent  components  of  a  waveform,  we  consider  here 
the  relative  difference  in  their  power  using  an  expression 
for  the  ratio  of  noncoherent  to  coherent  backscatter.  The 
expression  is  dependent  upon  the  radar  system  parameters 
as  well  as  the  radar  height  above  the  surface  and  the  rough 
surface  statistics.  At  nadir,  the  expression  is 

^  =2(A/)-/3^f  T  -ML  (20) 

where  k  =  2ii/A,  ko-A  and  /3pif  is  the  one-way  half  beam 

width  (/3^pif  is  equal  to  cx/H).  Note  that  the  ratio  becomes 
independent  of  the  surface  reflection  coefficient  when  this 
approximation  is  used. 

Figure  7-8  puts  into  perspective  the  relative  contribu¬ 
tions  of  coherent  and  noncoherent  backscatter  to  the  maxi¬ 
mum  power  in  a  peaked  signal.  The  upper  panel  shows  the 
ratio  for  a  satellite-borne  altimeter,  and  the  lower  for  an 
airborne  altimeter,  calculated  from  Equation  (20).  In  each 
case,  the  ratio  is  plotted  against  A:  <7,  with  isolines  of  Af.  The 
figure  indicates  that  either  component  of  backscatter  can 
dominate  the  signal  peak,  depending  on  the  roughness  of 
the  surface  and  the  altitude  and  system  point  target  re¬ 
sponse  of  the  radar.  For  the  regime  Act<  1 .0  and  kl  <  7.0,  the 
scattering  should  be  predominantly  coherent  according  to 
Kirchhoff  theory.  Such  surfaces  have  been  measured  for 
artificial  ice,  labeled  (a),  and  naturally  occurring  sea  ice  in 
the  Gulf  of  Bothnia,  labeled  (b).  Such  surfaces,  however,  are 
rarely  encountered  in  dynamic  regions,  where  typical  rms 
height  and  correlation  lengths  exceed  0.3  and  7.0  cm, 
respectively  (i.e.,  ko  >  1.0  and  kl  >  7.0).  Based  on  field 
measurements,  the  noncoherent  component  approaches 
and  can  exceed  the  coherent  component  in  the  roughness 
regimes  of  the  Beaufort  and  Labrador  Seas,  labeled  (c)  and 
(d). 
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In  Figure  7-8(a.i,  the  relevant  .sy.sten)  pararneiei are 
selected  for  the  EKS-1  altimeter.  Re,sults  indicate  (hal  it  i.s 
likely  that  a  dominant  coherent  component  will  be  ob,-<erved 
in  the  waveform  peak  over  mo.sl  con.soliduted  ice  .suriaees 
except  those  such  as  in  rough  ice  regime  (di  In  reginu  s  la  u 
(b I,  and  (Cl  the  coherent  component  i.s  more  than  •'>  dli  above 
the  noncoherent  component.  To  contrast  with  the  .satellite 


ka 

Fig.  7-8.  Incoherenticoherent  ratio  for  (ai  spacebome  (ER8-1 
altimeter) and  (b)  airborne  (RAL  altimeter)  platforms.  The  labeled 
areas  correspond  to  ice  roughne.ss  regimes  measured  in  .several 
experiments:  1  =  CRRFLEX,  0.29  <  ka  <  O..08.  2.9  <  M  <  8.7; 
2  =  BEPERS,  0.5b  <  ka  <  0.97,  \.2  <  kl  <  19.5;  ,3  Beaufurt'OO. 
0.85  <ka<  1.^  6,5.8  <kl  <  34.7;  and  4  =  HMEX’87, 1.19<^(7<  128, 
11.6  <  kl  <  66.5. 
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situation,  Figure  7-8(bi  demonstrates  that  the  simple  in¬ 
crease  in  the  antenna  i’LF  beam  width  of  the  airborne 
altimeter  results  in  increasing  dominance  of  difluse  or 
noncoherent  backscatter.  Kvidently,  only  in  the  smoother 
regimes  of  ice  roughness  will  the  coherent  component  ap¬ 
proach  or  exceed  the  noncoherent  returns.  This  may  occur 
over  calm  ocean  surfaces  in  marginal  ice  zones,  over  frazil 
slicks  or  smooth  new  ice  sheets  such  as  nilas,  or  even  during 
summer  melt,  when  maximum  melt  pond  coverage  is 
achieved.  Under  calm  wind  conditions  melt  ponds  provide 
extremely  effective  specular  reflectors. 

Typical  waveforms  from  the  RAL  altimeter  show  dis¬ 
tinct  examples  where  the  peaked  component  in  the  wave¬ 
form  does  dominate  the  diffuse  component  iDrinkwatcr, 
1 987 ;  1 991 1,  even  in  circumstances  where  the  majority  of  ice 
is  relatively  rough.  This  cannot  easily  be  explained  by  the 
simplified  arguments  oresented  above,  because  the  ice 
surface  in  those  cases  is  assumed  to  be  uniform  in  its 
characteristics.  Obviously  there  are  situations  whore  phase- 
coherent  returns  can  occur  from  rough  surfaces  with  mixed 
reflection  coefficients  and  roughness  scales.  This  idea  is 
applied  to  a  case  of  quasispecular  RAL  altimeter  signals 
where  the  ratio  of  coherent  to  noncoherent  powers  is  of  the 
order  of  25  dB  [Drinkw'ater,  1987;  19911.  Such  signals  are 
shown  to  originate  in  areas  of  smooth  open  water  between 
floes  (such  areas  have  a  reflection  coefficient  about  10  dB 
greater  than  that  of  ice),  and  areas  of  new  ice  growth.  For 
mixtures  of  specular  ocean  and  diffuse  ice  scattering,  the 
areal  jiroportion  of  smooth  ocean  required  to  produce  this 
type  of  signal  is  less  than  1  %  1  Drinkwater,  1991 1.  At  least 
in  the  marginal  ice  zone,  the  dominant  scattering  mecha¬ 
nism  therefore  appears  to  be  coherent  reflection  from  sn.ooth 
water  between  floes,  or  from  new  ice  growth  where  centime¬ 
ter-scale  waves  are  damped  out.  Similar  examples  of 
extremely  strong  signals  exceeding  30  dB  (and  often  reach¬ 
ing  40  dB  or  more)  have  been  observed  in  satellite  altimeter 
measurements  by  both  Ulander  |1987b;  1991 1  and  Laxon 
(1989).  Robin etal.  il983]calculatedthatonly0.01%ofthe 
surface  within  the  PLF  need  be  smooth  to  produce  a  char¬ 
acteristically  peaked  satellite  waveform. 

Given  the  importance  of  open  water  to  waveform  shape 
and  power,  it  is  necessary  to  understand  under  what  condi¬ 
tions  areas  of  open  water  might  produce  coherent  returns. 
Long  wavelength  swell  is  almost  completely  damped  out 
within  a  few  tens  of  km  of  the  ice  edge.  Shorter  wind  waves 
are  damped  in  a  short  distance,  on  the  order  of  km,  but  may 
be  regenerated  by  wind  in  areas  of  open  water  [Wadhams  et 
al.,  19881.  Wadhams  11983]  calculated  the  short  wave¬ 
length  roughness  generated  within  the  pack  using 


where  U,  is  the  friction  velocity  determined  from  wind 
speed,  X  is  fetch,  and  g  is  the  acceleration  due  to  gravity. 
The  standard  deviation  of  surface  height  versus  fetch  for 
various  wind  .speeds,  calculated  from  liquation  (21),  is 


shown  in  Figure  7-9.  Figure  7-9  illustrale.s  how  wind  .speed 
can  determine  whether  an  area  of  open  water  ol  a  given 
width  will  appear  coherent  to  the  altimeter  For  a  fre¬ 
quency  of  13  GHz,  this  require.s  that  O' be  le.ss  than  2.88  mm 
to  meet  the  Rayleigh  criterion,  or  less  than  0.72  mm  to  meet 
the  more  stringent  Fraunhofer  criterion.  Although  tfe.se 
relationships  are  onlv  approximate,  they  at  least  give  a 
first-order  indication  of  the  circumstances  under  which 
open  water  within  the  pack  is  likely  to  produce  a  near- 
specular  return. 

To  conclude,  the  scattering  mechanism  for  the  narrow- 
peaked  waveforms  obseiwed  by  altimeters  has  been  a  rather 
controversial  issue  in  the  past.  There  are  not  enough 
measurements  of  the  surface  roughness  of  ice  at  mm  scales 
to  provide  a  solid  basis  for  theoiy.  Investigators  have 
arrived  at  different  conclusions:  that  the  waveforms  origi¬ 
nate  from  true  specular  reflection  (Shapiro  and  Yapiee, 
1975;  Eom  and  Boerner,  1986]  or  from  quasispecular 
noncoherent  scattering  IRobin  et  al.,  1983;  Rapley,  1984; 
D  ’nkwater,  1987|.  The  latter  explanation  assumes  that 
tiie  rms  surface  slope  is  .small  (loss  than  the  angle  sub¬ 
tended  by  the  PLF),  such  that  the  angular  extent  of  the 
scattering  is  confined  to  within  the  altimeter  angular  re¬ 
sponse,  It  turns  out  that  the  effective  back.scattering 
coefficient  for  the  two  cases  in  fact  agree,  as  noted  by  Robin 
et  al.  [1983 1.  This  implies  that  it  is  impossible  to  unambigu¬ 
ously  determine  the  scattering  mechanism  from  waveform 
data  alone  [Ulander  and  Carlstrom.1991 1. 

One  approach  to  resolving  the  quertion  for  sea  ice  is  to 
compile  probability  density  functions  of  echo  energy  ovei 
example  areas.  If  the  scattering  is  noncoherent,  the  prob¬ 
ability  density  functions  (pdfs)  will  result  in  a  Rayleigh 
distribution.  If  indeed  coherent  returns  are  occurring,  the 
pdfs  will  be  more  Gaussian  in  nature.  The  compilation  of 
such  pdfs  should  be  undertaken  for  a  variety  of  instru¬ 
ments  of  differing  operating  parameters  (frequencies  and 
beam  widths)  in  order  to  establish  the  principal  effects. 
Alternatively,  it  should  be  possible  with  an  airborne  -system 


Fig.  7-9.  Standard  deviation  of  surface  height,  ,ersus  fetch  for  wind 
speed.s  from  1  to  1 6  m/s.  The  cutoffs  for  the  Rayleigh  and  the  more 
stringent  Fraunhofer  criteria  for  a  smooth  .surface  are  shown. 
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to  overfly  the  same  surface  at  different  altitudes.  This  will 
enaole  the  range  dependency  in  the  scattering  signature  to 
be  evaluated  and  thus  the  scattering  mechanism  to  be 
resolved  completely.  Comparing  airborne  and  helicopter- 
borne  results  is  an  alternative  method  of  analysis. 

Note  that  the  dependence  on  altitude  in  Equations  (11) 
and  (!2i  implies  that  it  may  be  difficult  to  compare 
spaceborne  and  airborne  measurements  if  a  coherent  com¬ 
ponent  is  present.  Typically,  the  altitude  of  an  airborne 
altimeter  is  two  orders  of  magnitude  less  than  that  of  a 
satellite,  which  reduces  the  coherent  component  from  that 
measured  by  a  satellite  altimeter  by  about  20  dB,  The 
noncoherent  component,  on  the  other  hand,  is  unchanged. 

Realistically,  in  considering  coherent  reflection  one  must 
acknowledge  that  n’any  returns  with  different  tracker 
biases  or  tracker  jitters  are  added  together  to  produce  a 
mean  satellite  waveform.  This  inevitably  results  in  a 
waveform  which  is  the  equivalent  of  the  random 
superposition  of  mixed  coherent  and  noncoherent  compo¬ 
nents  of  the  return  signals.  The  fading  of  signals  occurring 
between  independent  along-track  views  of  Ihe  surface  fac¬ 
ets  will  inherently  lead  to  some  mis-superposition  of  pulses 
and  thus  noncoherent  averaging  of  the  signals  during  the 
tracking  process. 

7.4  Thkoky  Vkk.si'S  Mka.sukement 
7.4.  J  Surface  Meciaurements 

To  improve  the  retrieval  of  geophysical  information  from 
altimetry,  it  is  necessary  to  better  understand  what  an 
altimeter  is  sensing.  Answers  to  questions  such  as — Where 
is  the  signal  originating?  and  What  are  the  effects  of  snow 
depth  and  ice  type  on  backscatter?  can  be  provided  by 
surface  measurements.  Near-surface  backscatter  mea¬ 
surements  from  the  Fram  Strait  marginal  ice  zone  (MIZ) 
have  been  acquired  at  frequencies  of  5.25,  9.6, 13.6,  and 
16.6  GHz  using  a  frequency-modulated  continuous-wave 
radar  scatterometer  with  pencil-beam  antenna  patterns 
(two-way  beam  widths  of  1 .5°  to  4.5°).  The  radar  scattering 
coefficients  were  determined  by  comparing  surface  returns 
with  returns  from  calibration  targets  of  known  radar  cross 
section.  These  observations  were  made  from  a  helicopter  in 
conjunction  with  surface  observations  that  allow  backscat¬ 
ter  to  be  related  to  a  particular  ice  type  or  to  ocean  condi¬ 
tions  [Onstott,  1991 ;  Onstott  et  al . ,  1 9871.  Also  given  in  this 
section  is  a  synopsis  of  recent  measurements  in  the  Weddell 
Sea. 

7.4,1, 1.  Vertical  incidence  scatterometer  data  from  the 
Fram  Strait  marginal  ice  zone.  The  four  major  ice  types 
pre.sent  when  the  mea.surcments  were  made  in  midsummer 
were  multiyear  (MY ),  thick  first-year  (TFY),  medium  first- 
year  (MFY),  and  thin  first-year  (ThFY).  Snow  depths  for 
MY  and  FY  ice  by  midsummer  were  about  half  of  what  they 
were  at  the  beginning  of  summer.  This  is  important  in  tuat 
thick  ice,  such  as  TFY  and  MY  ice,  will  have  enough  snow 
cover  to  mask  much  of  the  ice  surface  from  microwave 


radiation  for  al  least  half  of  the  iwo-raonlh-long  .suninier. 
Data  were  actjuired  liy  tlying  scans  from  a  pu.sition  over  the 
ocean  about  20  km  from  lb<‘  ice  edge,  acros.s  the  ice  edge,  and 
into  the  MIZ  until  large  floe.s  were  encountered  iaii  addi 
tional  20  km).  One  such  Iran.sect  is  .shown  in  Figure  7-10, 
Illustrated  in  the.se  data  is  the  transition  from  the  strong 
backscatter  of  the  open  ocean.  t(.  a  region  of  calm  water  (  5 
km  in  e.vtcnt '  immediately  adjacent  to  the  ice  edge,  through 
a  well-delineated  and  compact  ice  edge,  into  a  region  ,viih 
bunds  of  similarly  sized  ice  floes  and  varying  ice  concentra¬ 
tion,  and  then  into  a  region  ol'OO'i  ice  concentration  and 
'argc  multiyear  floes. 

Using  the  open  ocean  signature  as  a  reference,  the  area 
of  dead  vvater  near  the  ice  edge  produced  an  enhancement 
of  about  6  dB-  Backscatter  from  ice  is  dramatically  smaller 
(by  8  to  21  dBt.  in  contrast  to  the  usual  ob.servations  by 
airborne  and  satellite  altimeters  moving  from  ocean  to  MIZ. 
This  difference  may  be  explained  by  the  much  larger  Fresnel 
zone  of  airborne  and  satellite-borne  instruments,  which  can 
encompass  flat  areas,  with  a  resulting  dominance  of  coher¬ 
ent  reflection.  As  the  ice  concentration  increa.ses  deeper  in 
the  MIZ,  a  lower  average  return  is  observed.  The  strongest 
returns  are  associated  with  the  water  between  floes.  These 
returns  are  specular  in  nature  and  are  characterized  by  a 
cross  section  at  vertical  of  greater  than  25  dB  and  a 
reduction  in  intensity  of  almost  five  orders  of  magnitude  by 
25°.  These  results  support  the  assertion  that  peaked 
returns  in  airborne  and  satellite  data  come  from  smooth 
water  between  floes.  The  range  of  scattering  coefficients  is 
also  of  interest.  The  range  for  the  open  ocean  is  about  6  dB. 
whereas  in  the  MIZ  it  is  about  40  dB.  The  high  range  in  the 
MIZ  is  due  to  the  presence  of  areas  of  open  water,  which 
produce  a  variety  of  microwave  signatures  influenced  in 
large  part  by  their  position  from  the  ice  edge  ( Figure  7-11). 
Outside  and  within  the  ice  edge,  variations  in  ocean  surface 
roughness  are  caused  by  complicated  interactions  between 
wind,  waves,  currents,  and  fresh  water  input.  These  inter¬ 
actions  are  driven  by  effects  related  to  cold  air  from  the  pack 
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Fig.  7-10,  Measurements  from  a  singleradar .scattcromotertransect 
acros.s  the  ice  edge  on  July  .5,  1984  at  13.6  (IHz  and  a  vertical 
incidence  angle. 
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crotioinj^  tro  il  the  celd  tieeaii  uce  and  sea  suri.ue  U‘nip4.-ra 
ture  oft)  Ui  2  5"  Ci  near  the  at*  I'dtjf  to  the  warnu'r  optni 
Ofeamseasurt'ace  temperature  ot  -I  t  ('  >,  One  result  ol'thts 
movement  ofcoki  air  is  the  area  ofealm  water  adjacent  to 
the  ice  edee. 

The  transition  from  the  open  ocean  to  the  calm  water 
within  the  MIZ  is  very  rapid  m  this  casi  (t  occurs  vvuhin 
tiOO  m,  and  it  produces  kiroe  level  shitLs  in  backscatier,  '1  his 
sharp  delineation  was  observed  often  during  M  IZKX’tsd  and 
’84.  Ambiguities  or  errors  in  detecting  the  ice  edge  with 
altimetry  may  occur  if  the  ocean  near  the  ice  edge  damped 
by  a  cold  air  mass  is  mistaken  as  an  ocean  surface  damped 
by  a  thin  layer  of  ice.  In  the  satellite  case,  the  ice  edge  is 
marked  by  a  similar  increa.se  in  variability  but  higher 
back.scatter  values  over  the  ice. 

The  MIZ  is  a  dynamic  region.  DilTu.se  .scattering  from 
deformed  ice  is  ob.served  to  produce  returns  12  dH  Ivelow 
that  of  flat  ice  Kelurns  from  water  in  the  shadows  of  floes 
aa  2  to  5  dli  above  typical  ocean  returns  from  within  the 
MIZ.  Backscalter  modulations  are  also  found  to  lx*  a 
function  of  ice  concentration  and.  to  a  lesser  extent,  fltx' 
size.  Based  on  other  ob.serv  atioiis  by  Onstotl,  the  backscat- 
ter  response  for  winter  is  anticipated  to  lx>  very  similar  The 
distribution  of  features  within  the  bounclarie.s  of  a  floe  may 
also  have  an  important  impact  on  the  backscalter  respon.se. 
Pressure  ridges,  melt  pond.- ,  and  the  blocks  laat  pile  around 
the  floe  edges  require  con.iideralion.  New  pre.ssure  ndges 
may  form  at  floe  boundaries  and  are  difl'erentiated  from 
older  pres.sure  ridges  that  may  be  covered  with  a  layer  of 
wcl.snow.  P'dges  result  in  weak  back.scatter,  with  the  more 
recently  formed  ridges  producing  the  more  difTu,sc  scatter 
response.  Pool.s,  often  small  and  protected  from  wind  by- 
undulating  ice,  were  found  to  produce  backscalter  cross 
.sections  of  a  magnitude  .similar  to  that  of  shadowed  water 
signatures.  On  July  5.  when  these  observations  were  made. 

20  nm  From 
Ice  Edge 

10  nm  From 
Ice  Edge 

0.5  nm  From 
fee  Edge 


1  nm  in  MIZ 


10  nm  in  MIZ 


In  Pack  Ice 


Floe  Shadow 
Effect 

10~o'“  15.0  2ab  25.0"  30.0’”"’ 35^ 

Radar  Scattering  Cross  Section.  dB 

Fig.  7-1 1 .  Liquid  oce.in  radar  scattering  cross  .sections  shown  a.s  a 
function  of  position  from  the  ire  i-dge  at  1  3.f>  <}li/.  and  a  vertical 
incidence  angle. 
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mimnial  floe  are.i  wa.s  covered  l>\  ineh  poiids  Me!*,  poud 
formaluiii  w.i.s.  isuwi-ver.  under  r.ipid  ;u  celt  t .ilion  le,.!  ir, 

■  fuly  tl.  I'livi-r.ige  on  .MV  floe.-,  wu.s  e.-jtjtu.iti  d  .n  .li.out  . 

7  -1  /  2.  let  lypf  (Hid  feature  lii-enii.-.ruituai  t'ru’e 
Strati  ttifunurernfnts.  During  sutfilio;!  ,  Uif  leilet '.)\ it >  .0. 
nadir  of  F^'  and  M^'  ice  is.  deternnneci  f,y  the  pnn  im-  o; 
moist  snow,  the  areal  extent  of  nu-h  pond-,  -md 
.Melt  poiuls  inc.rea.se  itie  retleclivily  of  a  floe  .ind  topogra 
phy  acts  to  reduce  the  intensity  of  the  reflection  At 
13  6  till/.,  the  refleclivilK-s  fur  the  three  pninarv  ne  !vi.es 
are  witlun  a  very  narrow  i  .'idBrange  fabieV  i  .sfiow  >  ttiat 
the  backscatter  re.'-ponse  follow.-,  in  Lirgi  part,  -urfaie 
topography.  Thin  fir.-t  year  ice  i-  cli.tr. leten.-'ed  fiv  tfo 
ihmriesl  .snow  cover  inuan  of  4  cm  and  e  taljuiar.  vvitii 
little  if  any  deformation  it  produces  the  largest  r<.  flei !  iv  ilv 
of  the  three  ice  lyjx-s  I’hvsicallv,  .51F'V  ice  j.s  sm.ii.it  to 
ThFY  ice.  except  that  it  has  a  thu  kei  snow  cover  -  no  .m  of 
10  cmi,  MFY  and  MY  .share  a  sitniiar  sign.itur<  .  of  .m 
ice  sheet  with  a  very  thick  snow  layer  Thi-  d..ita  at  !  .3  ti  t'.H/ 
for  snow-covored  .MY  apjx-ar  high  at  fir.'t  giant  o,  hut  upon 
further  examination  they  are  consi.stenl  witti  the  oth.er 
oh.servations  .shown  in  Table  7  !  A  .MF\'  flix  produced  a 
signature  -3  dB  lower  than  that  of  ThFY.  w  here  .snow  vva-- 
about  1 0  cm  On  the  same  .MFY  fltx*.  however,  ifien  was  .i 
regio’  where  snow  had  accumulated  to  14  cm  I  lore  tin- 
reflectivity  increa.sed  to  8  if  dB.  Thi.s  is  identical  to  the 
reflectivity  ob.served  on  the  MY  ice  and  within  1 ,2dB  of  t  hat 
of  ThF'Y  ice.  Thick  snow  on  tiie  .MY  ice.  rather  than  the 
snow-ice  inlerface.  dominates  the  MY  .signature 

While  backscatter  for  winter  and  summer  FT  ice  i.  very 
similar,  significant  difTerences  in  the  .seasonal  respon.se  of 
MY  ice  are  observed  The.se  are  rel.ited  to  the  moisture 
contained  m  the  snow  and  ic<  sheet  [On.slotl  et  al  .  !9S7i 
Moisture  limits  the  penetration  into  the  ice  .sheet,  therefiv 
reducing  or  negating  the  volume  seattcTing  that  dominates 
the  backscatier  respon.se  of  MY  ice  dunng  winter  .Xt 
vertical  incidence,  this  cau.ses  a  9-dB  signature  enhance¬ 
ment  Uil  13.6  GHzi  for  MY  ice  dunng  summer  when  com¬ 
pared  to  the  winter  response.  Measurcmcnl.s  of  snow 
moisture  show  that  the  bulk  wetness  of  the  snow  pack  is 
about  5  5'7  by  volume,  and  the  majority  of  the  .snow  pack ’s 

TABLE  7-1  Badar  scattering  ccK-fncient,-  at  vertical  inriticncc 
angles expres.sed  in  decilxds.  Measurement.-  were  nhtained  .hiiv  ii. 
1984.  in  the  marginal  ice  zone  lOnstott,  !  991  !. 
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old  and  has  a  donsily  of  400  to  500  kg  m  *  oi  al  . 

1987  i.  Thf  penetration  depths  ealculated  for  iht^se  condi¬ 
tions  are  about  4.5  cm,  1 .8  cm,  1 .5  cm,  and  1  .:5  cn.  fur  5.25, 
9.6, 13.6,  and  16.6  GHz,  respectively  The  .snow  pack  on  .\IY 
ice  up  to  midsummer  may  be  40  cm.  Note  that  al  1 .3.6  GHz, 
only  1  to  3  cm  of  .snow  rs  required  before  the  .snow  layer 
completely  dominate.s  the  back.scatier  re.sponse  since  for 
each  penetration  depth  the  incident  signal  experience.s  a 
9-dB  round-trip  loss. 

7. 4. 1.3.  Shiphorne  Ku  barul  radar  uhscnatiotr-t  tn  Ihr 
Weddell  Sea.  Ku-band  radar  data  were  acquired  during  a 
transect  of  the  K.V.  Polarstern  through  the  Weddell  Sea 
during  the  austral  winter  of  1989 !  Lytle  et  al.,  1992i,  A  VV 
polarized  network  analyzer-ba.sed  radar  deveU)[H“d  at  the 
University  of  Kansas  made  ob.se  rvations  at  incidence  angles 
spanning  O'  to  303  Most  radar  data  were  acquired  on  ice 
floes  against  which  the  ship  was  moored.  .Measuremenl.s  of 
physical  properties  of  the  sea  ice  included  surface  rough¬ 
ness,  snow  cover  thickness,  snow  and  icc  density,  and  gram 
size  and  .salinity.  Ice  freelxiard  data  were  also  collected  by 
drilling  numerous  holes  through  the  ice  and  measuring 
w'ater  level  relative  to  the  ice  surface. 

Flooding  at  the  ice-snow  interface  .seems  to  be  a  charac 
teristic  phenomenon  of  the  western  Weddell  Sea  i  Lange 
and  Eicken,  1991;  Wadhams  et  al..  198"!,  This  process  is 
apparently  associated  with  enhanced  snow  accumulation 
depressing  the  ice  surface  below  .sea  level.  Flooding  changes 
density  and  salinity  in  the  lower  parts  of  the  snow  layer, 
which  in  turn  affects  the  radar  response.  Backscatter  is 
plotted  against  freeboard  and  parameterized  against  slush 
layer  presence  in  Figure  7-12.  Backscatter  i,s  elevated  when 
a  slush  layer  is  present,  presumably  because  of  the  en¬ 
hanced  dielectric  contrast.  Flooding,  then,  may  be  another 
mechanism  for  situating  a  relatively  bright  target  within 
the  radar  footprint.  The  possibility  of  flooding  must  be 
taken  into  account  when  attempting  to  interpret  altimetry 
in  terms  of  ice  concentration  or  type.  Laboratory  measure¬ 
ment  to  further  study  the  importance  of  flooding  are  de¬ 
scribed  in  Chapter  9. 
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Fig.  7-12.  Average  ice  freeboard  versus  a  parameterized  against 
the  presence  of  a  slush  layer. 
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Vanou.-.  work  .MciiitMi-  ci  .li  ,  i ‘U'li  !  li  iuk  w ,it<  i  F*-'. 
Fedor  and  Wal.-.h,  !  ;  .h  i.->  dut  iiHi.  lil.-d  lii.-  t  tsar  .u  li  r  i,- 

lics  <if  airhuriif  altnin-u-r  i  i  iii>x!,un  ■  u.i 

ic«'  and  -water  Si-ierai  ui  tin  »•  .•ludji  -  i'laweii  i  l  ,ii  i  5 
Cowan  and  .Squire,  !‘,)->4  iii><  unu  nl  re-  ull:-  in-ni  .s  pi,- 
launch  c.'  ipaigii  tor  KK.8-1  ,  m  whu  l:  tlie  lOl.  .dt ma  lt  i 
acquired  dal.i  m  expenineiit.a  de'.iM'ii  t...  mda.ite  l‘a 
charac ten.stics of. -vca  ice rt-Ic. .int  Ui.i!' inii-Ii.  i  puf-t  .ni.tly  se- 
.Airtiorne  aktme!  a!  lierierally  ii.e.  e  v.  uier  iM  ani-"  Uku;  ~:t\ 
ellile  altimeters,  iii  order  to  reduce  tlie  iinp.ict  o!  .UileiUi.! 
oO'-pointliig  cau.seil  liy  aircr.ift  inotion  .A  w  ide  he.un  !«  .tij  - 
to  different  waveform  char.ic'enstic-  ituin  .ire  i4i,--<  r-. a.  d 
with  a  ---atellite-tiorne  m.slruint:rit .  a.s  wa.-  di.'tu,'-.ei:i  in 
Section  7.3  2.  The  Mirlace  area  .-.enia  d  bv  .m  .iirbtirrit- 
iruslrument  i,s  n.'-ualiy  a  facior  of  ii.ni  sm.iih  i  llian  tor  a 
satellite  altimeter,  while  it  .--ubti  ntj.-.  a  wider  riinge  at 
incidence  angk-.s.  i'his  i.-  a  cntu.i!  diirerence.  bet  ui,'*-  ui 
area.s  consisting  of  nu-Xture.-'  of  ice  and  w.iier  tin  re  i-. 
inherently  a  larger  .irea  of  water  .ivailabk-  to  the  .-.iieHue 
altimeter  for  .sc.ittenng  from  nadir  tlian  i.--  availalile  (n  ihe 
lower  altitude  airlKirne  instriirnenl  .Another  eonsi-qut net 
of  the  airlwrne  altimeter.-  wider  lie.nn  is  that  bat  kstitUer 
fall-olT  with  inciiierue  angle  can  be  measured 

7 .4.2. 1.  IVfij  efarm  example.^.  Figure  7- 1  3  is  an  exatnfde 
of  the  operation  of  the  liAI-  aUimeter  over  a  varii.-ty  rif  sea 
ice  conditions  in  tlie  marginal  ice  zone  of  ihe  Frau.  ,81,  .iii 
Con-sccutive  waveforms  exhibit  dramatic  charuje,-  in  pro¬ 
portion  as  the  instrument  pa.s.scs  over  large  ice  floes  and 
area.sofopon  water  in  the  pack  Note  that  .some  waveforms 
.show  .saturation.  One  mleresling  finding  illustrated  h\  the 
waveform  labeled  a  is  that  as  the  altimeter  pa.-se.s  over  a 
lead  only  a  few  meters  wide,  the  altimeter  responds  ex¬ 
tremely  sharply  if  the  lead  is  at  nadir.  Fedor  and  Walsh 
119881  noted  the  .same  phenomenon  in  data  from  the  Ik  au- 
fort  Sea.  Since  the  power  reflection  coefTK  ient  of  water  is 
generally  1 1  dB  higher  than  that  for  snow  -covered  sea  ice, 
the  water  between  ice  floes  dominate.s  the  hackscattered 
signal  near  nadir,  Olhiu  ohseiwalion.s  from  thi.-  dam  set 
show  that  when  the  sea  ice  fills  the  footprint,  waveform 
shape  becomes  suggestive  of  ocean  returns  as  the  signal 
level  is  modulated  by  the  type  and  .surface  roughnc.ss  of  sea 
icc  present.  F'cdor  et  al.  (1 9891  note  typical  peak  backscal- 
ter  strengths  of  around  30  to  37  d  B  over  areas  of  new  .sea  ice 
and  around  1 0  dB  over  deformed  F'A'  ice  in  the  Beaufort  Sea. 
These  observations  show-  close  correspondence  with  the 
findings  of  Drinkwatcr  11991  !  over  F’Y  ice  in  the  FTam 
Strait.  In  the  latterstudy,  accompanyingsurface  measure¬ 
ments  [Tucker  et  al.,  1987]  led  to  attempts  to  da.ssifv  icc 
types  using  the  altimeter  waveforms.  The  study  revealed 
four  distinct  surface  clas.scs  based  on  surface  roughness  ( or 
fall-off  of  back.scalter  with  incidence  angle i  and  normal 
incidence  back.scatter  coefTicient.  Tho.se  corre.sponded  in 
near-coincident  SAR  images  and  aerial  photographs  with 
smooth, bare  F’Y icc  with  minimal  ridging;  snow-covered  F'Y 
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Fig,  7-13.  MIZEX'84  RAL  airborne  altimetry  over  mixed  ice 
concentrations  in  the  Fram  Strait.  Waveforms  show  received 
power  in  range  gate  bins  IDrinkwater,  1.9871 


ice;  rotted,  rough  conglomerate  ice  floes;  and  MY  ice  with 
large-.scale  surface  roughness. 

In  ice-water  mixtures,  the  rapidity  with  which  the 
power  decays  after  the  peak  is  largely  a  function  of  the  wave 
slope  f  see  Section  7.3.4),  Some  areas,  such  a.s  for  waveform 
6,  provide  non-fetch-limited  regions  of  almost  open  water 
within  the  pack.  Wind-generated  gravity  and  capillary 
waves  develop  rapidly  in  the  absence  of  floes  (unlike  ate). 
The  waveform  trailing  edge  responds  dramalieally  to  the 
presence  of  these  waves  and  particularly  to  the  higher  wave 
.slopes.  The  waveforms  become  more  ocean-like,  albeit  with 
a  relatively  high  normal  incidence  backscatter  coefticient 
due  to  the  lack  of  large  gravity  wave.s. 

Volume  scattering  cannot  be  neglected  when  old  ice  has 
low-density,  bubbly  upper  layers  (Chapter  2).  Volume 
scattering  leads  to  a  backscatter  re.sponse  function  that 
tapers  much  more  gradually  with  incidence  angle.  It  is 
sugge.sted  as  a  mechanism  for  some  Fram  Strait  observa¬ 
tions  where  neither  the  rm.s  surface  roughness  nor  the 
residual  aircraft  motion  can  account  for  the  spreadingof  the 
waveform  leading  edge.  Together  with  the  large-.scalc 
surface  roughness  observed  on  vast  low-salinity  MY  does, 
volume  scattering  appears  to  result  in  rms  tracker  delay 
variability  or  height  fluctuations  of  between  0  and  4  m 
IDrinkwater,  1991 1.  The  snow  depth  on  these  floes  aver¬ 
aged  28  cm,  and  ridge  .sail  heights  varied  between  3  and 


5  m.  Ill  contra.'-l.  '.ivi-i  .('Y  lin  nu  .i.-ui'.-'i 

-surfaev  height  van.ibshu  w.i.--  K;.---.  i;i,'.!i  1  .nn-u-! ,  atid  .-.ah 
tracking  i  xcursiois-  were  not  noted  For  the  ^.3ni•  et'i  >  t.i 
occur  111  satellite  data,  the  radar  would  !ik<  iv  h.c.e  li, 
directly  over  a  v.rst  .\1Y  ice  iloi-  or  iter'err  Tin  -itu.ilion 
would  be  eiicuuiitered  in  the  .\i  ctic  diu  in;.;  mien  oral; tu. a  e- 
and  in  the  Antarctic,  vvfiere  va,-.t  tabular  m  ixu;-'-  an  edieo 
ob.served  by  allimelry  I.McliUvrt'  ami  ('iidlip  Itts;"  ui 
1989;  and  llaw  km.s  et  al  ,  1991 

7-/2  2.  Ire  runrrrtt  ra!  !>n:  i,''t  r  iiu'io'i  j ’',trr  i  a  [  r!<  r o.'i 
scn  atinns.  .'.ttempt.--  have  been  made  to  u>e  .liliiiu  tn  lor 
ice  concentr.itioii  nu-.isureineiil:-  Tlie  iiiodels  den  rii«-d  m 
Section  7.3  3  yield  fractioMal  area.'-  olAiiiuii!  li  ice  ami  ii 
water  facets,  if  reOectnin  coefticient  .iiid  sna-  are  -  piei  slu  t* 
These  do  not.  however.  iiKorpor.ite  tile  .^lope  disli  ibo' ion, 
which  1.S  portrayed  for  both  ice  floe-  and  ilie  water  U  iwi .  n 
in  the  heuristic  model  in  Figure  7- 1  i.  ilu:-  modfi  likistrate.' 
the  decline  in  signjfR-.inl  wave  height  cau.-ed  liv  ,-i  .t  hi 
wave  allenuation  i  Raplev,  1  9.s.| ;  and  the  a:-,-.oi  i.ited  dei  line 
in  wave  slope  a.-  short  wave  length  wave.-  .ire  fireierenlianv 
damped.  The  resulting  modulation  of  tfie  w.iveforni  ))i,>wer 
directly  from  ice  concenir.itioii  will  therefore  wiry  .i-  .i 
function  ofboth  rro.-  wave  height  and  fr.tctional  coverage  of 
.sea  ice 'based  on  the  premi.se  that  the  returned  power  in  tiie 
waveform  decays  a.s  the  ratio  of  low-rellectivity  n  e  to  high- 
reflectivity  water  increasesi  Figure  7-1  4  sliows  ideah/cii 
waveforms  at  concenlration.s  of  ().  20.  40,  60,  so.  and  !  00' . 
for  a  given  wave  regime  In  each  ca.se  the  tola!  height 
distribution  is  a  comlnned  effect  of  wave  height .  ice  .surlao 
roughness,  and  penetration,  w  hde  the  formof  t  fir  waveform 
trailing  edge.s  is  a  result  ed  tlie  wave  .-lope.-  ohsi.’rved 
between  ice  (Iocs. 

Using  altimeter  data  ove  r  ocean  suiTares  inside  t  he  M IZ 
edge,  Dnnkwater  1 1  9871  shows  that  the  angular  hackseat- 
ter  .scattering  function  i.s  predicted  hv  noncoherent  geo¬ 
metric-optics  theeiry.  The  mode!  fit  to  the  data  give,-  an  rms 


Fig.  7-1 4.  Wnvefoniis  generated  by  heuristir  model  orwrive  (leray 
with  attenuation,  and  families  of  w  a  vefomi.sas.socia  ted  with  varying 
ICO  concentration:  0,  20,  40,  60.  80,  and  lOO'i  inrinkwater.  1!)H7! 
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slope  12.8-'.  This  is  consistent  with  waves  measured  in  the 
ice  edge  for  a  fairly  well-developed  sea.  With  increasing 
distance  into  the  MIZ,  observations  show  a  corresponding 
decrease  of  rms  slope  to  around  3  ’  at  40  km  from  the  ice 
edge.  The  progressive  damping  of  ocean  waves  as  they  pass 
into  the  MIZ  results  in  more  and  more  spiked  echoes  from 
the  water  between  ice  ftoes,  as  shown  in  F'igure  7-14.  Thus, 
waves  and  the  wave  spectrum  in  the  MIZ  mutually  interact 
w  ith  the  shape  and  strength  of  altimeter  signals  and  must 
be  considered  when  attempting  to  derive  concentration 
from  waveforms.  In  an  analysis  ofthe  relationships  between 
backscattered  power  and  ice  concentration  in  the  Greenland 
Sea,  Drinkwater  IT 991 1  investigates  airborne  Ku-band  sig¬ 
nals  for  which  there  are  accompanying  SAR  and  photogra¬ 
phy  data.  To  obviate  the  requirement  for  compensating 
wave  height  information,  altimeter  waveforms  are  chosen 
beyond  the  point  at  which  they  are  modulated  by  waves. 
These  are  subsequently  compared  writh  digitized,  simulta¬ 
neous  aerial  photographs  over  a  track  length  of  130  km. 
Both  the  peak  backscattered  power  and  integrated  power  in 
the  trailing  edge  of  the  waveform  are  shown  to  be  related  to 
the  ice  concentration.  However,  because  of  the  saturation 
effects  found  in  these  data,  the  peak  power  could  not  be  used 
for  accurate  ice  concentration  retrievals.  The  waveform 
integral,  on  the  other  hand,  declines  monotonically  with 
increasing  ice  concentration  and  is  reasonably  robust.  Re¬ 
gression  analysis  shows  that  the  relationship  between 
these  two  variables  had  a  correlation  coefficient  of  -0.81, 
and  65'/c  ofthe  signal  energy  variance  is  accounted  for  using 
this  measure  alone.  The  rms  residual  error  is  13%,  and  the 
remainder  of  the  variance  is  accounted  for  by  varying 
surface  conditions  due  to  local  wind  roughening,  varying  ice 
types,  ice  roughness,  and  other  factors. 

7.4.3  Satellite  Measurements 

A  common  problem  in  the  analysis  of  altimetry  is  to  find 
high-quality  surface  truth  data  with  which  to  determine  ice 
conditions  during  the  sensor  overflight.  For  corroborating 
the  interpretation  of  satellite  altimetry  the  most  successful 
data  source  has  been  SAR  [Swift  et  al.,  1985,  Ulander 
1987a,  Fettereret  al.,  1991,  Ulander  and  Carlstrbm,  1991  j. 
Examples  of  investigations  that  rely  on  SAR  are  given  in 
this  section.  Comparisons  of  altimetry  with  imagery  from 
passive  sensors  is  reserved  for  Section  7.4.4. 

7.4.3. 1.  Echo  waveform  simulations  compared  with 
measurements.  We  noted  earlier  that  it  is  impossible  to 
di  ^ermine  unambiguously  the  scattering  mechanism  from 
a  single  narrow-peaked  waveform.  However,  a  series  of 
such  waveforms  may  provide  this  information  as  the  altim¬ 
eter  moves  from  a  smooth  area  into  an  adjacent  area  with 
diffusely  .scattering  characteristics.  For  true  specular  re¬ 
flection,  the  returned  power  essentially  originates  from  the 
first  Fresnel  zone.  This  implies  a  horizontal  resolution 
according  to  Equation  (9),  which  is  almost  an  order  of 
magnitude  finer  than  the  pulse-limited  resolution  for  an 
ERS-1  class  altimeter.  Hence,  the  echo  waveform  will,  in 


the  specular  case,  respond  more  quickly  to  a  change  ifoni  a 
narrow-peaked  to  a  diffuse  type  of  scattering. 

This  observation  has  been  used  to  determine  the  scatter¬ 
ing  mechanism  of  narrow-peaked  waveforms  over  sea  ice  in 
the  Baltic  Sea  1  Ulander  and  Carlstrdrn,  1991 1.  A  section  of 
Geosat  altimeter  data  was  extracted  over  a  newly  frozen 
lead  surrounded  by  very  deformed  ice.  The  former  produced 
narrow-peaked  waveforms  and  the  latter  essentially  dif¬ 
fuse  waveforms.  A  pulse  echo  .scrie.s  wa.s  then  .simulated  by- 
integrating  over  the  illuminated  surface  al  each  platform 
position  along-track,  and  an  echo  waveform  was  produced 
by  averaging  over  one  hundred  pulse  echoes  The  surface 
composition  is  shown  in  Figure  7-1 5(a).  It  was  derived  from 
a  high-resolution  airborne  SAR  image  that  delineated  the 
two  surface  types.  Figure  7-1 5(b)  shows  the  Geosat  normal 
incidence  backscatter  coefficient  and  two  simulations  that 
assume  a  coherent  and  noncoherent  mechanism,  respec¬ 
tively.  In  the  noncoherent  simulation,  the  young  ice  was 
represented  by  a  model  according  to  Equation  (1 7),  whereas 
in  the  coherent  simulation  all  the  returned  power  origi¬ 
nated  from  the  first  Fresnel  zone  with  uniform  weighing. 
Note  the  dose  agreement  between  the  Geosat  data  and  the 
noncoherent  simulati,..!.  The  model  parameters  in  Equa¬ 
tion  (17)  were  chosen  as  F  =  3.3%  and  =  0.068*  (for  Z>  - 
5.5  m).  The  coherent  simulation,  on  the  other  hand,  shows 
a  much  higher  sensitivity  to  the  changing  sui  !ace  charac¬ 
teristics  than  does  the  Geosat  data.  This  particular  ex¬ 
ample  shows  that  the  narrow-peaked  waveforms  over  the 
frozen  lead,  corresponding  to  a  high  backscatter  coefficient 
of  32  dB,  were  indeed  noncoherent  in  nature.  The  back¬ 
scatter  value  over  the  lead  in  Figure  7-1 5(b)  is  somewhat 
less  due  to  distortion  by  the  pulse  echo  simulation  and 
waveform  averaging. 

The  remaining  issue  to  resolve  is  whether  the  noncoherent 
scattering  is  in  accordance  with  the  geometrical  optics  or 
flat  patch  model.  Both  models  predict  a  backscatter  coeffi¬ 
cient  of  the  same  functional  form  and  have  essentially  the 
same  effective  backsca tiering  coefficient  for  narrow-f>eaked 
waveforms.  Note,  however,  that  the  cause  for  the  angular 
response  is  quite  different.  In  the  geometrical  optics  model 
it  is  determined  by  the  rms  surface  slope,  since  the  reflect¬ 
ing  specular  points  are  infinitely  large  compared  to  the 
wavelength;  in  the  flat  patch  model  it  is  determined  by  the 
flat  patch  diameter  since  the  patches  are  assumed  horizon¬ 
tal  (i.e.,  with  zero  surface  slope).  There  are,  as  yet.  no  field 
measurements  which  could  help  resolve  this  question.  In 
general  it  would  require  that  the  joint  height  probability- 
density  function  be  determined  for  scale  lengths  between  1 
cm  and  100  m.  Hence,  our  effort  in  the  future  must  turn 
toward  more  sophisticated  field  measurements  of  the  sur¬ 
face  roughness. 

7.4.3.2.  Ice  type  signatures  in  satellite  measurements. 
With  an  approach  quite  different  from  that  outlined  in 
Section  7.4.2  for  airborne  altimetry,  Chase  and  Holyer 
(1990]  extend  the  measurement  of  ice  type  and  concentra¬ 
tion  to  satellite  altimetry.  Their  work  uses  a  linear  un  mixing 
model,  suggested  by  the  fact  that  waveforms  from  the  MIZ 
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Fig.  7-15,  <a)  Simuiated  surface  composition  derived  from  an 
airborne  SAR  image  of  Baltic  Sea  ice,  and  (b)  simulation  result 
compared  with  Geosat  measurements  [Ulander  and  Carlstrbm, 
1991!. 

and  pack  are  usually  the  result  of  backscatter  from  mix¬ 
tures  of  ice  and  water  in  the  footprint.  Each  waveform  is 
represented  as  a  multidimensional  vector  that  can  be  ex¬ 
pressed  as  a  linear  combination  of  end  member  vectors. 
Analysis  consists  of  identifying  the  end  members  and  asso- 
ciatingeach  with  a  surface  type.  Although  the  model  can  be 
adapted  for  nonlinear  combinations  of  waveforms,  Chase 
and  Holyer  argue  that  it  is  not  nece.ssary  for  the  data  they 
examined.  Using  manually  interpreted  airborne  passive 
microwave  imagery  collected  in  the  Greenland  Sea  within 
hours  of  a  Geosat  overpass,  Chase  and  Holyer  found  that 
the  analysis  of  the  concentration  of  three  combinations 
of  types  (open  water  and  brash  ice,  FY  and  old  ice,  and 
young  and  new  ice)  agreed  with  that  from  the  passive 
microwave  analysis  to  within  4.6%  over  the  experiment 
area.  Interestingly,  the  shape  of  the  end  member  paradigm 


vcmvefonn.-i  .sugge.st.s  that  concentration  in  thi.-.  anaic.-  i.--  i.- 
directly  pruporlional  to  peak  power  and  wavel'orin  nitiyir  li. 
Possible  explanation.s  iiir  tln.-^  contra.sl  with  the  airhorne 
case  are  the  dilfcrencc.s  in  .scale  and  range  di.-ciia.-ed  in 
Section  7.3.2. 

Alost  investigations  oi  the  ability  of.-atellite  alt  nneti.  r 
measureinenl.s  to  distingui.sh  between  broad  ice  thu  kne.-s 
cla.s.ses  and  to  detect  w  ater  o|)eniiig.i,  ,-turting  w  ith  that  of 
Dwyer  and  Godin  1 1 9S0 1.  have  only  usetl  the  altimeter  data 
in  a  relative  sense.  This  i.s  because  tlie  data  lack  ahsohitc 
radiometric  calibration.  It  should  be  ernphusi/ed  that  it  i> 
only  pos.sihle  to  ohtiun  an  elTective  backscatter  coeiriciiait 
w'hen  the  angular  waveform  is  narrow  peaked.  By  a.s.sum- 
ing  a  particular  form  of  the  angular  resfronse.  it  may  then 
be  pos.sihle  to  obtain  a  be-tter  estimate  of  the  hack.scalter 
coefficient  at  normal  incideiiee.  This  inhimently  assiirni  s 
that  the  .sy.stom  impul.se  respon.se  i.s  suOiciently  narrow  to 
resolve  the  angular  scattering  .-esponse  and  that  the  scat¬ 
tering  mechanism  is  known. 

The  correlation  between  radar  hack.scalter  at  nonn.il 
incidence  and  Arctic  and  Baltic  winter  ice  types  has  been 
described  by  .several  inve.stigators  I  Ulander.  198Tu. 

Fedor  et  al.,  1989;  Ulander  and  Carl.stnim,  1991  f  lli.'-lo- 
grams  of  the  normal  incidence  backscatter  coenicienl  for 
some  of  these  ice  types  are  shown  in  Figure  7-16.  The 
normal  incidence  back.scatter  coefl'icient  decrea.ses  in  gen¬ 
eral  for  increasing  ice  thickness,  .such  that  new  ice  has  the 
highest  values  (30  to  40  dB)  and  .MY  ice  the  lowe.st  (10  to 
20  dB).  Very  deformed  ice  also  has  a  low  back.scatter 
coefficientdOdB).  Note  the  high  mean  values  for  al!  classes 
when  compared  with  the  surface  measurements  of  Sect  ion 
7.4.1.  These  results  suggest  that  broad  ice  clas.ses  can  he 
inferred  from  the  altimeter  backscatter  data,  provided  the 
surface  conditions  do  not  change  rapidly  along  track.  When 
the  latter  is  not  the  case  (for  instance  in  the  .MlZs  the 
signature  shows  a  large  variability,  w'hich  introduces  dif¬ 
ficulties  in  the  interpretation. 

Figure  7-17  shows  curves  of  the  backscatter  coeflicient 
and  the  echo  waveform  width  (a  measure  of  diffuse  scatter¬ 
ing)  for  a  Seasat  altimeter  track  in  early  October  off  the 
coast  of  Alaska  lUlander,  1987bl.  The  curves  are  overlaid 
on  a  Seasat  SAK  image  that  w'as  acquired  a  few  hours 
earlier.  Narrow,  peaked  waveforms  dominate  over  the 
pack,  except  where  the  nadir  track  crosses  tw  o  large  floes. 
Figure  7-1 8  shows  a  similar  comparison  ofGeosat  altimetry 
and  airborne  SAR  in  the  Greenland  Sea  IFettcrer  el  al , 
1991).  Here,  however,  the  high-frequency  variability  in 
backscatter  is  much  greater.  Backscatter  is  from  the 
0.1-second  waveforms.  Note  that  AGC  i.s  a  .somewhat 
smoothed  version  of  peak  power.  The  voltage  proportional 
to  attitude,  or  VAIT  parameter,  is  a  measure  of  power  in 
later  gates.  Low  values  of  V'ATT  (plotted  every  .<econd ) 
indicate  peaked  returns.  The  parameters  in  Figure  7-18(0 
delineate  zones: 

1 )  The  MIZ,  where  backscatter  rise.s  from  9  to  -30  dB  as 
the  icc  edge  iscrossed.  Waveforms  are  strong  and  peaked, 
suggesting  reflection  from  w’ater  between  the  .small  lloo-s, 
which  give  a  bright  return  in  the  SAR  image,  Back.scatter 
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Fig.  7-16.  Backscatter  histograms  for  (a)  Beaufort  Sea  ice  from 
Seaaat  d-ata  [Ulander,  1988],  and  (b)  Baltic  Sea  ice  from  Geosat 
data  [Ulander  and  Carlstrom,  ©  1991  IEEE]. 
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drops  slightly  over  an  ex.pan.-r.o  o!  open  watt  r  .A-vera!  kilo 
meters  wide,  where  wind  may  roughen  the  .^urt'aee. 

2,1  The  pack,  where  vanabilily  i.s  high,  Baekstatter  ir 
somewhat  low  er  in  the  we.slern  half,  where  there  os  a  larger 
percentageof  rough  ice.  Back.scalier  dip.s  over  two  large  old 
floes,  which  appear  quite  bright  in  the  S.-VR  image 

3]The  fast  ice,  where  back.scatter  reache.s  it-s  lowest  level 
and  wavefonns  are  not  peaked  but  are  more  characteristic 
of  (X’ean  returns.  Thi.s  i.s  evident  in  the  waveforims  plotted 
in  Figure  7-18(  a  >. 

While  changes  in  the  nature  and  variability  of  wave¬ 
forms  are  revealing,  caution  should  lx‘  taken  in  interpreting 
the  power  and  shape  of  wa  veform.s  from  an  inhomogeneous 
area  such  as  this  in  terms  of  ice  roughne.ss  or  type.  Instru¬ 
mental  effects  and  rapid  changes  in  surface  characlerist  ie.s 
(described  in  Section  7.2,5)  drastically  alter  individual 
averaged  waveforms. 

7.4.4  Satellite  Altimetry  Compared  With  Infrared.  Vinihle, 
and  Passive  Microwave  Imagery 

In  Sections  7.4.2  and  7.4.3  aerial  photography  or  SAR 
imagery  was  used  to  gain  insight  into  scattering  mecha¬ 
nisms  and  to  test  the  possibility  of  deriving  concentration 
and  type  estimates  from  altimetry  on  a  case-by-case  basis. 
Further  applications  arise  from  the  synergistic  use  of  radar 
altimetry  with  imaging  sensors.  For  example,  if  a  certain 
ice  type  can  be  identified  in  imagery,  this  may  aid  greatly  in 
the  interpretation  of  the  radar  altimeter  data,  which  may 
then  be  used  to  measure  surface  elevation  and  roughness 
(derived  from  waveform  leading  edge  width).  This  section 
compares  altimetry  with  wide-swath  imagery  on  the  regional 
and  global  scales  required  of  operational  algorithms.  Such 
comparisons  enhance  our  understanding  of  the  operation  of 
both  sensors  while  aiding  the  development  of  altimeter 
data  algorithms.  Immediate  applications  lie  in  the  area  of 
algorithm  validation.  For  example,  the  altimeter  has  shown 
great  promise  in  validating  ice  extent  mapping  by  passive 
microwave  instruments.  In  addition,  ERS-1  altimeter  sea 
ice  products  can  be  validated  using  data  from  the  ERS-1 
Along  Track  Scanning  Radiometer. 

7.4.4. 1.  Waveform  parameterization.  Figure  "-1 9  shows 
the  two  parameters  computed  from  altimeter  waveforms  to 
quantify  the  response  over  sea  ice  in  the  following  analysis. 
The  first,  SIGPK,  represents  the  peak  backscatter  mea¬ 
sured  in  the  return  echo.  The  second,  SIGTD,  is  computed 
as  the  difference  between  the  average  power  in  eight  early 
and  eight  late  range  gates  in  the  altimeter  return.  SIGTD 
is  essentially  a  measure  of  the  degree  of  peakiness  in  the 
return  echo. 

7.4.4.2.  Observations  of  a  compact  ice  edge  in  the 
Greenland  Sea.  The  Advanced  Very  High-Resolution  Radi¬ 
ometer  ( AVHRR)  onboard  the  NO  AA  series  of  pol  ar  orbiting 
satellites  provides  visible  and  infrared  imagery  with  a 
swath  width  of 2900  km  and  nadir  resolution  of  1 .1  km.  It 
is  routinely  used  for  operational  ice  analysis.  Figure  7 -20(a) 
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Fig.  7-17.  Seasat  altimeter  measurements  in  the  Beaufort  Sea  compared  with  Seasat  SAR  for  October  3, 1978  [Ulander,  1987b]. 


shows  four  Geosat  altimeter  tracks  overlaid  on  a  charmel  2 
(visible  band)  image  of  the  Greenland  Sea.  Both  altimeter 
data  and  imagery  were  acquired  on  May  10, 1987.  Param¬ 
eters  SIGPK  and  SIGTD  are  derived  for  every  waveform. 
There  is  therefore  a  value  every  0.1  seconds,  or  every  0.7  km 
along  the  ground  track.  A  sharp  rise  in  both  parameters  is 
observed  close  to  the  ice  edge  visible  in  the  imagery.  In  some 
cases,  notably  tracks  B  and  C,  a  rise  in  SIGPK  is  observed 
outside  the  ice  edge  and  before  any  noticeable  rise  in 
SIGTD.  The  cause  of  this  is  unknown,  but  the  examples 
demonstrate  the  sensitivity  of  the  altimeter  to  the  presence 
of  sea  ice.  A  similar  rise  in  backscatter,  related  to  the 
damping  of  waves  due  to  the  proximity  of  sea  ice,  is  seen  in 
the  surface  measurements  of  Section  7.4.1.  The  highest 
strength  (with  values  up  to  40  dB)  and  most  peaked  echoes 
are  observed  just  inside  the  ice  edge,  which  is  normally 
populated  by  small  (<1 00  m)  ice  floes.  Further  into  the  pack 
the  waveform  strength  and  peakiness  decrease  and  become 
more  variable.  This  is  attributed  to  a  less  heterogeneous 
surface  caused  by  the  presence  of  larger  ice  floes.  In  some 
areas,  such  as  where  track  A  crosses  track  D,  cloud  cover 
can  confuse  identification  of  the  ice  edge.  In  such  cases  the 
altimeter  may  be  used  to  better  identify  the  loca  .ion  of  the 
ice  edge. 


7. 4.4.3.  Observations  of  a  diffuse  ice  edge  in  the  Greenland 
Sea.  Figure  7-20(b)  shows  three  altimeter  tracks  overlaid 
on  an  AVHRR  channel  4  (infrared)  image  acquired  on 
March  18, 1987.  In  this  instance  an  off-ice  wind  results  in 
a  more  diffuse  ice  edge  than  is  observed  in  the  May  10 
image.  The  two  altimeter  tracks  labelled  A  and  B  show  a 
gradual  increase  in  echo  strength  over  some  distance  with 
little  increase  in  echo  peakiness.  This  is  attributed  to  a 
more  gradual  damping  of  ocean  waves  than  is  observed  in 
the  May  10  image.  Again  the  strongest  and  most  peaked 
returns  are  observed  further  into  the  ice  pack. 

Of  particular  interest  in  this  image  are  the  two  large 
floes  observed  close  to  the  coast  and  transected  by  altimeter 
track  B.  The  floes  are  clearly  delineated  by  low  strength, 
ocean-like  echoes.  Unlike  other  areas  of  the  pack,  where 
smooth  ice  or  open  water  can  occur  within  the  altimeter 
footprint,  over  these  vast  ice  floes  the  altimeter  is  observing 
only  a  rough  ice  surface.  Such  a  phenomenon  may  make  it 
difficult  to  discriminate  between  the  return  from  open 
water  and  that  from  a  large  floe.  However,  because  such 
surfaces  fill  the  footprint  for  several  waveform  integration 
times,  they  do  provide  the  possibility  of  measuring  meter- 
scale  roughness  and  even  ice  freeboard.  The  AVHRR  image 
was  acquired  some  13  hours  after  altimeter  pass  B.  The 
southward  offset  of  the  two  ice  floes  with  respect  to  the 
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Fig.  7-18.  (a)  Geosat  altimeter  waveforms  from  April  13, 1987,  with  a  (b)  near-coincident  (within  1  hour)  airborne  SAR  image  of  the  East 
Greenland  Sea  at  72°  N.  White  lines  on  the  image  mark  the  approximate  position  of  the  9-km-wide  range-window-limited  footprint,  (c ) 
Measures  of  peak  power,  a’(0),  average  power,  AGC,  and  waveform  shape,  VATT,  are  plotted,  Along-track  registration  of  SAR  image 
and  altimeter  data  is  estimated  to  be  within  2  km  [Fetterer  et  al.,  1991 ). 


signal  in  the  altimetry  can  be  attributed  to  motion  in  the 
East  Greenland  current. 

A  similar  ocean-like  signal  is  observed  over  the  fast  ice, 
with  a  sharp  decrease  in  echo  strength  at  the  fast  ice 
boundary  followed  by  a  slight  step  increase  in  echo  strength 
further  into  the  coast.  Again  this  is  due  to  the  altimeter 
observing  only  a  homogeneous  rough  ice  surface.  The 
delineat'on  of  two  zones  within  the  fast  ice  suggests  that  the 
altimeter  is  discriminating  between  an  inner  zone  of  fast  ice 
that  has  survived  the  winter  intact,  and  an  outer,  rougher 
zone  that  has  been  broken  up  at  some  time  by  a  winter 
storm.  (Similar  zonation  in  the  returns  from  fast  ice  is 
apparent  in  Figure  7-18). 


7. 4.4.4.  Midwinter  image  of  the  Kara  Sea.  Figure  7-20(c) 
shows  three  altimeter  tracks  overlaid  on  an  AVHRR  chan¬ 
nel  4  image  of  the  Kara  Sea  acquired  on  February  21, 1986. 
Waveform  echo  data  are  not  available  for  this  period  and 
therefore  the  AGC  values  are  displayed.  In  this  thermal 
image,  darker  areas  indicate  higher  surface  temperatures 
and  therefore  thinner  ice.  To  the  right  of  the  island  Novaya 
Zemlya  an  area  of  relatively  new  ice  has  formed.  On  the  far 
right-hand  side  of  the  image  a  further  large  frozen  lead, 
some  20  km  across,  is  observed  close  to  the  continental  land 
mass.  For  this  image,  echo  strength  shows  a  strong  corre¬ 
lation  with  ice  age.  New,  smooth  ice  produces  stronger 
returns  than  older,  rougher  ice.  This  phenomenon  is  ex¬ 
actly  what  would  be  expected  when  considering  normal 
incidence  radar  backscattcr  over  surfaces  with  varying 
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Fig.  7-19.  Wavefonn  parameters  SIGPK and  SIGTD  are  measures 
of  the  peak  power  in  a  waveform  and  of  the  peakiness,  or  rate  of 
backscatter  fall-off  with  incidence  angle,  in  the  waveform.  SIGTD 
is  calculated  by  taking  the  difference  in  average  power  between 
eight  early  and  eight  late  range  gates. 

degrees  of  roughness.  The  highest  echo  strengths  are 
observed  over  the  narrow  lead,  which  is  slightly  darker 
than  the  ice  just  east  of  Novaya  Zemlya.  Within  the  older 
ice  some  cracks  are  observed  to  correlate  with  individual 
peaks  in  the  altimeter  AGC  value.  In  some  cases  peaks  are 
observed  in  AGC  that  do  not  correlate  with  any  cracks 
visible  in  the  imagery.  Smooth  ice  or  open  water  not 
resolved  in  AVHRE  imagery’  are  believed  responsible  for 
these  peaks.  These  observations  demonstrate  the  sensitiv¬ 
ity  of  altimeter  data  to  areas  of  thin  ice  and  leads.  The 
mapping  of  leads  has  great  importance  to  sea  ice  research, 
as  leads  modify  ocean-atmosphere  heat  transfer  in  the  high 
Arctic. 

7.4. 4. B.  Midwinter  image  of  the  Chukchi  Sea. 
Figure  7-20(d)  shows  three  Eiltimeter  tracks  overlaid  on  an 
AVHRR  channel  4  image  of  the  Chukchi  Sea  acquired  on 
March  14, 1988.  All  three  tracks  show  numerous  peaks  in 
the  echo  strength,  especially  near  coasts  and  islands.  Here, 
leads  are  most  likely  to  occur  when  winds  blow  older  ice 
away  from  the  coast,  edlowing  new  ice  to  form.  Judging  by 
high  SIGPK  and  SIGTD  values ,  a  large  area  of  new  ice  may 
be  north  of  Wrangel  Island,  although  this  is  difficult  to 
observe  in  the  imagery.  This  again  shows  the  usefulness  of 
altimetry  as  an  adjunct  to  visible  or  infrared  imagery. 

7.4.4.6.  Altimetry  coincident  with  passive  microwave 
imagery.  Satellite  passive  microwave  data  has  been  widely 
used  for  mapping  global  sea  ice  cover.  Satellite  altimetry 
provides  the  only  practical  alternative  means  for  global 
synoptic  monitoring  of  sea  ice  extent.  The  merging  and 
comparison  of  data  from  these  two  sensors  is  therefore  an 
obvious  means  to  improve  understanding  of  both  [Laxon 
and  Askne,  1991]. 

In  the  following  analysis  we  examine  data  from  a  merged 
product  produced  using  data  from  the  Geosat  altimeter  and 
the  SSM/I  flown  on  board  the  Defense  Meteorological  Sat¬ 
ellite  Program  (DMSP)  satellite  series.  The  five  brightness 


temperatures  recorded  by  the  SSM/I  for  the  day  and  loca¬ 
tion  of  each  altimeter  data  point  are  appended  to  the  Geosat 
Geophysical  Data  Record  (GDR).  The  spatial  coincidence  is 
therefore  within  25  km  (the  gridded  resolution  of  daily 
average  SSM/I  brightness  temperatures)  and  the  temporal 
coincidence  within  one  day.  The  merged  data  set  can  then 
be  used  to  directly  compare  the  ice  extent  mapped  by  the 
two  sensors.  The  SSM/I  data  are  processed  using  the  NASA 
Team  algorithm.  The  algorithm  used  for  detecting  the  ice 
edge  in  the  Geosat  data  record  uses  the  standard  deviation 
of  the  10  surface  height  values  recorded  in  the  GDR  over 
1  second.  Over  sea  ice,  the  disruption  of  the  tracker  (de¬ 
scribed  in  Section  7.2)  results  in  a  significant  increase  in  the 
standard  deviation  of  the  height  measurement  [Laxon, 
1990].  Other  algorithms  have  used  the  rise  in  AGC  and  the 
drop  in  SWH  to  detect  the  ice  edge  and  the  penetration  of 
swell  in  the  MIZ  [Rapley,  1984]  and  the  rise  in  the  ratio  of 
AGC  to  VATT  that  occurs  as  the  edge  is  crossed  [Hawkins 
and  Lybanon,  1989]. 

Figure  7-21  shows  a  series  of  altimeter  tracks  for 
August  3, 1987  overlaid  on  a  polarization  difference  image 
generated  from  the  85  GHz  channel.  Although  the  85  GHz 
channel  suffers  from  some  weather  mtamination,  the  ice 
cover  is  reasonably  visible.  A  data  j .  mt  is  plotted  in  white 
where  both  sensors  indicate  that  ice  is  present.  Where  only 
the  altimeter  indicates  ice  the  data  point  is  plotted  in  red. 
Where  only  the  SSM/I  indicates  sea  ice  is  present  the  data 
point  is  plotted  in  green.  If  neither  sensor  indicates  that  ice 
is  present,  a  blue  data  point  is  plotted.  Both  sensors  show 
an  ocean  signal  in  most  areas  away  from  the  ice  boundary, 
with  the  exception  ofsome  green  points  north  of  Scandinavia 
where  the  SSM/I  incorrectly  identifies  some  areas  as  being 
ice  covered.  It  is  well  known  that  the  SSM/I  maps  some 
areas  of  open  water  as  ice  due  to  weather  effects.  In  the 
Beaufort  and  Chukchi  Seas,  various  areas  ofboth  green  and 
red  points  are  observed.  The  red  points  appear  to  lie  near 
the  ice  boundary,  while  the  green  points  occur  within  the 
pack  (in  one  case  over  what  might  be  an  area  of  open  water). 
Small  areas  of  red  points  are  also  observed  in  the  Kara  and 
Greenland  Seas.  These  may  delineate  ice  cover  below  the 
lowest  concentration  detectable  using  SSM/I. 

In  order  to  look  at  the  Geosat  altimeter  and  SSM/I 
differences  on  a  larger  scale,  two  17-day  merged  data 
sets  were  compiled,  one  during  midwinter  and  one 
during  midsummer  for  both  hemispheres.  The  data  plots 
produced  from  these  data  sets  do  not  show  any  points  where 
both  instruments  indicate  no  ice  to  be  present.  This  is  to 
permit  clearer  delineation  of  areas  of  disagreements.  Fig¬ 
ure  7-22(a)  shows  data  for  August  1987  in  the  Northern 
Hemisphere.  The  amount  of  ice  cover  below  72°  N  (the 
northern  limi  t  of  coverage)  is  limited,  but  significant  clus¬ 
ters  of  red  points  are  observed  in  the  Greenland,  Beaufort, 
Barents,  and  Kara  Seas.  At  lower  latitudes,  such  as  around 
the  British  Isles,  a  large  number  of  green  data  points 
illustrate  the  problem  of  weather  effects  in  the  SSM/I  data. 
Figure  7 -22(b)  shows  the  corresponding  image  for  February 
1988.  During  maximum  extent  there  is  much  better  agree¬ 
ment  overall  between  the  altimeter  and  the  SSM/I.  The  few 


Fig.  7-20.  (a)  An  AVHRR  visible  image  of  the  Greenland  Sea,  with  altimeter  waveform  parameters  overlaid  along  the  Geosat  altimeter 
track  (Scoresby  Sound  is  near  the  center  of  the  image);  (b)  as  in  (a),  but  with  a  less  compact  ice  edge;  (c)  infrared  band  image  of  the  Kara 
Sea,  with  Geosat  AGC  values  plotted  every  0.1  s  along  the  nadir  track;  and  (d)  infrared  band  image  of  the  Chukchi  Sea,  with  waveform 
parameters. 


exceptions  are  in  the  Chukchi  Sea,  around  165°  E,  and  in 
the  Greenland  Sea,  where  some  green  points  are  observed. 
These  points  may  be  attributable  to  areas  of  fast  ice  and 
vast  floes  which,  as  seen  in  the  comparisons  with  AVHRR 
data,  result  in  an  altimeter  signal  that  is  difficult  to  distin¬ 
guish  from  open  water. 

Figure  7-22fc)  shows  the  Southern  Hemisphere  during 
February  1988.  There  is  a  clear  geographic  zonation  of 


gre  n  and  red  points.  Since  the  altimeter  samples  the  ice 
edge  at  random  points  during  the  17-day  period,  the  fact 
that  the  areas  of  disagreement  are  geographically  corre¬ 
lated  strongly  suggests  some  geophysical  cause.  Red  points 
in  particular  are  observed  mainly  in  the  Weddell  and  Ross 
Seeis,  with  green  points  occurring  in  the  Bellingshausen 
Sea  and  along  the  eastern  coast  of  the  Antarctic  continent. 
The  cause  of  this  discrepancy  is  not  understood,  but  the 
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Fig.  7-21.  An  August  3, 1987,  SSM/I  85-GHe  polarization  differ¬ 
ence  image. 


implications  for  ice  edge  mapping  are  significant,  since  the 
altimeter  is  routinely  identifying  an  ice  edge  nearly  1 00  km 
outside  that  identified  by  SSM/I  in  some  areas.  The  green 
points  on  the  eastern  pa-t  of  the  continent  might  be  attrib¬ 
utable  to  fast  ice,  but  those  in  the  Bellingshausen  Sea  are 
more  difficult  to  interpret.  In  contrast  to  the  Northern 
Hemisphere  weather  effects  do  not  seem  to  result  in  so 
many  green  points  over  the  open  ocean.  Figure  7-22(d) 
shows  the  corresponding  image  for  August  1987.  Again 
there  is  better  overall  agreement  during  the  winter  period, 
with  very  few  red  points  occurring  anywhere.  Some  areas 
of  green  points  are  observed  mainly  in  the  Bellingshausen 
Sea.  The  reason  for  this  is  not  fully  understood.  It  may  be 
attributable  to  areas  of  nearly  100%  ice  cover,  resulting  in 
an  ocean-like  altimeter  signal. 

7,5  Summary 

Potential  users  of  radar  altimetry  will  find  sources  for 
data  from  past,  present,  and  future  missions  listed  in 
Massom  [1 991  ] ,  with  a  treatment  of  the  applicability  of  the 
data  to  polar  studies.  Methods  for  extracting  various  ice 
parameters  from  the  data  record  have  been  suggested  in 
this  chapter.  Additional  parameters  and  alternative  ways 
of  deriving  parameters  can  be  found  in  the  literature  (e.g., 
Rapley  etal.  119871). 

Generally,  ice  parameters  are  arrived  at  by  first 
retracking  and  editing  the  altimeter  data  record  and  then 
deriving  waveform  parameters  that  are  linked  to  the  de¬ 
sired  ice  parameters.  The  processing  required  to  derive  a 
waveform  parameter  and  the  strength  of  theory  linking  it  to 


an  ice  puranuder  vary  I'landiT.  iur  inslaitc.',  <  .dt  - 
backscatter  from  wavefonn.s  and  cla.-'.-.c-.s  lee  on  the  b.c-i-  oi 
bnckscatter  lUlander.  1991  i.  Drinkwatcr  |lvi91  1  rind.-  the 
theoretically  inver.-e  relation.-hip  between  the  intcitral  u!' 
waveform  power  and  ice  conceiilration  borne  out  in  inea 
surement.s.  Cha.seand  Ilolyer  119901  periorm  a  aojihi,-;!!' 
cated  analy.si.s  of  waveform  .shape,  and  eniiiincally  link 
shape  to  different  ice  type.s  and  concentraiions.  Gtln.r 
algorithms  use  parameter.s  that  are  more  dost  ant  !y 
derived  from  waveform  power  atid  shape.  For  instame, 
SWH  has  been  used  to  mca.surc  swell  penelrtition  of  ilie 
MIZ,  and  AGC  has  been  u.sed  to  mark  the  ice  I'dge  n-ee 
Section  7.-1. 4).  While  AGC  and  SW’li  can  be  prnce.ssi  <1 
quickly  andare  more  readily  avtiilable  than  the  waveforrtis 
themselves,  theyare  the  re.sult  of  averaging  wavelbrm.stitid 
therefore  have  disadvantages.  Yet  another  apjtroach  is  to 
link  variability  in  waveform  parameter.s  with  ico  condition.- 
Variability  in  waveform  power  aid.s  in  the  delineation  of 
zones  in  Figure  7-18,  while  the  standard  deviation  of  the 
height  measur'^ment  is  used  in  Sect  ion  7  1  4  to  locale  the  ice 
edge. 

With  the  exception  of  ice  edge,  it  ha.'-  not  been  demon¬ 
strated  that  any  ice  parameter  can  be  reliably  retrieved 
from  altimetry  on  more  than  a  ca.se  .study  ha.sis.  The 
research  reported  here  strongly  indicates  that  this  will 
change.  The  best  hope  for  progress  in  thi.-  direction  may  lie 
in  seeking  an  empirical  connection  between  waveform  pa¬ 
rameters  and  ice  conditions  usingglobal  satellite  data  sets. 
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